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11  AeSTRACT 


HILRuW,  the  second  United  States  underground  nuclear  explosion  on  Amchitka  Island,  was 
.•■n  order  of  magnitude  larger  than  LONG  SHOT,  the  first  test  there.  Magnitude-yield  scaling 
between  the  two  for  both  surface  waves  and  body  waves  followed  th'oretical  predictions  well 
and  agreed  with  Nevada  Test  Site  (NTS)  scaling.  A  computed  locatio.’'  of  HILROH  using  LONG  SHOT 
travel-time  anom.tlies  was  only  1.2  km  in  error  horizontally  and  11.8  km  vertically,  a  sub¬ 
stantial  improvement  over  a  location  using  no  anomalies.  The  presence  of  direct  shear  waves 
and  Love  waves  for  HILROH  can  be  attributed  to  causes  other  than  tectonic  strain  release,  such 
as  mode  conversion,  crack  formation,  and  structural  features  of  the  surrounding  medium.  The 
detection  of  Love  waves  from  LONG  SHOT  and  the  MILROH  cavity  collapse  were  important  results 
bearing  on  the  determination  of  the  shear-generating  mechanism.  All  the  common  identification 
criteria  were  applied  to  MILROH  and  its  collapse:  location  and  depth,  H  vs  m.  ,  spectral  ratios 
complexity,  shear-wave  excitation,  and  radiation  patterns.  HILROH  appealed  to“be  typical  of 
explosions  while,  in  contrast,  the  collapse  was  much  like  an  earthquake. 
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ABSTRACT 


MILROW,  the  second  United  States  underground  nuclear 
explosion  on  Amchitka  Island,  was  an  order  of  magnitude 
larger  than  LONG  SHOT,  the  first  test  there.  Magni tude -yi el d 
scaling  between  the  two  for  both  surface  waves  and  body 
waves  followed  theoretical  predictions  well,  but  there  are 
significant  intercept  changes  in  the  magnitude-yield 
relations  from  the  Nevada  Test  Site  to  Amchitka  Island. 

A  computed  location  of  MILROW  using  LONG  SHOT  travel-time 
anomalies  was  only  1.2  km  in  error  horizontally  and  11.8  km 
vertically,  a  substantial  improvement  over  a  location 
using  no  anomalies.  The  presence  of  direct  shear  waves  and 
Love  waves  for  MILROW  can  be  attributed  to  causes  other 
than  tectonic  strain  release,  such  as  mode  conversion, 
crack  formation,  and  structural  features  of  the  surround¬ 
ing  medium.  The  detection  of  Love  waves  from  LONG  SHOT  and 
the  MILROW  cavity  collapse  were  important  results  bearing 
on  the  determination  of  the  shear-generating  mechanism. 

All'  the  common  identification  criteria  were  applied  to 
MILROW  and  its  collapse:  location  and  depth,  M^  vs  m^^, 
spectral  ratios,  complexity,  shear-wave  excitation,  and 
radiat-nn  patterns.  MILROW  appeared  to  be  typical  of  explo¬ 
sions  wnile,  in  contrast,  the  collapse  was  much  like  an 
earthquai'e. 
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INTRODUCTION 

On  23  October  1965  the  first  underground  nuclear  explo¬ 
sion  on  Amchitka  Island  in  the  Aleutians,  LONG  SHOT,  was 
detonated.  It  engendered  a  voluminous  analysis  of  distant 
sei smol ogi cal  data  and  significantly  increased  our  knowledge 
pertaining  to  location  and  identification  of  underground 
nuclear  tests,  as  well  as  our  general  knowledge  of  the  earth. 

On  2  October  1969  a  second  explosion,  MILROW,  was 
detonated  on  Amchitka  Island  with  a  yield  about  an  order  of 
magnitude  larger  than  LONG  SHOT.  The  purpose  of  this  report 
is  to  document  and  analyze  the  teleseismic  data  from  MILROW, 
mainly  that  recorded  by  the  LRSM  mobile  vans  deployed  at  ^’he 
time  and  by  the  three  V£LA  arrays,  UBO,  TFO,  and  l.ASA.  Some 
additional  data  taken  from  the  WWSS  network  of  the  NOS 
is  presented,  but  is  not  to  be  considered  complete  or  final 
at  this  time.  Important  studies  of  near-source  behavior 
are  being  conducted  by  government  agencies  and  private 
contractors;  their  results  have  not  reached  publication 
although  some  preliminary  reports  have  been  presented 
orally  (McKeown  et  al.,  1970;  Blackford  et  al.,  1970). 
Eventually,  we  hope  that  the  near-source  studies  can  be 
integrated  with  the  teleseismic  data  to  explain  the  MILROW 
source  dynamics  in  detail. 

Since  most  of  the  data  available  for  this  report  at  the 
Seismic  Data  Laboratory  and  through  the  WWSS  was  recorded  at 
sites  identical  to  those  for  LONG  SHOT,  since  MILROW  was 
detonated  only  2.3  km  from  LONG  SHOT,  and  since  the  overall 
quantity  of  data  is  reduced,  we  will  avoid  presentation  of 
any  analysis  which  would  merely  be  redundant  in  relation  to 
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the  vast  amount  pf  literature  previously  prepared  about  LONG 

\ 

SHOT.  In  particular,  little  new  sei smologi cal  information 
is  to  be  expected  from  the  P  and  PcP  signals  of  MILROW  because 
of  the  excellent  and  numerous  global  recordings  of  these  \ 
signals  from  LONG  SHOT. 

This  report  will  place  attention  on  the  presence  of 
shear  phases  which  were  not  seen  on  LONG  SHOT  recorr  ngs. 

The  magnitude-yield  scaling  \over  the  limited  range  ,f  80  to ^ 
1000  kt  will  be  ascertained  and  related  to  theoretical  pre-  ' 
(dictions.  The  occurrence  of  a  cons i deriabi e  collapse  of  the 

,  I 

■MILROW  cavity  two  days  afterward  produced  data  relevant  to 
the  determination  of  whether  MILROW  caused  any  tectonic  | 

istrain  release.  A'  thorough  application  to  MILROW  of  fcurrent 
discrimination  criteria  will  be  made  and  the  results  will 
be  compared  with  those  of  LONG  SHOT. 

Basic  spatial  and  temporal  location  |nformation  on 
MILROW  and  Hs  collapse  are  given  in  Table  I.  The  SDL's 
own  report  on  LONG  SHOT  (Lambert  et  al.,  1969)  is  the 

basic  reference  for  this  report,  and  it  is  hoped  that 

\ 

the  reader  will  ha\^|e  access  to  that  reference.  To 
eliminate  repetition,  we  will  usi^ally  omit  further  ref¬ 
erence  to  Lambert  et  al.;  and  data,  resuTts,  or  conclu¬ 
sions  pertaining  to  LONG  SHOT  discussed  in  this  report,  ' 

un^ss  attributed  to  another  source  or  newly  produced  in 
this  report,  are  taken  from  that  LONG  SHOT  repprt. 


\  ■ 
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THF  SOURCE  ENVIRONMENT 

I 

Physical  setting 

A  summary  of  gebl ogi cal  and  geophysical  data  concerning 

I 

the  Aleutian  arc  and  Amchit|<a  Island  in  particular  was  given 
in  the  LONG  SHOT  report  and  need  not  be  repeated  here.  One 
additional  noteworthy  reference  is  the  comprehensive  report 
onl  the  Aleutians  by  Anderson,  1970.  A  li*:ho‘10gic  log  at  the 
MILROW  site  (Snyder,  1969)  revealed  mainly  consolidated  vol¬ 
canic  breccia  beds  with  some  layers  of  andesite,  basalt,  and 
sandstone.  For  compari son ,  the  lithology  of  the  LONG  SHOT  site, 
2.3  km  away,  was  mostly  interbedded  tuff  and  breccia. 
Maqnitude-vs-yield 

LONG  SHOT  was  detonated  within  an  andesite  sill  while 
MILROW  was  detonated  in  breccia;  howevet,  at  these  deijjths  \ 
(2000-4000  feet)  the  overburden  pressure  and  degree  of  water- 
saturation  may  influence  the  medium  response  and  degree  of 
energy  radiation  more  than  the  exact  composition  of  the  rock. 
For  reference,  we  use  the  recent  magni tude-vs-yield  empi ri cal 
documentation  by  Evernden  (1970)  on  "hard  rock  and  mesa 
tuff"  events,  between  which  there  is  no  cl  ear, separati on . 

LONG  jSHOT,  reported  to  be  80  ikt  (USAEC,  1967)  1  has  an  (5.85) 
which  is  greater  by  about  three-tenths  of  a  unit  than  predicted 
by  events  detonated  in  these  media;  and  MILROW,  reported  to 
be  1  Mt  (USAEC,  1969),  also  has  an  mj^  (6.42)  greatei^  by  the 
same  amount  than  predicted  on  the  basis  of  these  other  events. 
All  of  the  events  for  hard  rock  and  mesa  tuff  in  Evernden's 
illustrations  are  from  the  NTS;  thus  the  effect  of  a 
combination  of  a  c|ifferent  source  environment,  different 
travel  paths,  and  different  recording  networks  on 
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the  magni tude-vs-yiel d  curve  is  apparently  to  increase 
body-wave  magnitude  a  few  tenths  of  a  unit  for  a  given 
yield  on  Amchitka  relative  to  NTS.  Orphal  et  al.  (1970) 
reported  that  the  close-in  response  at  1.0  cps  was  equal 
to  that  predicted  for  a  1  Mt  event  on  the  basis  of  NTS 
experience.  We  must  state  that  both  the  reported  yields 
for  LONG  SHOT  and  MILROW  are  pre-shot  figures,  subject 
to  any  errors  inherent  in  predicting  the  yield  of  a 
nuclear  device. 

Near-source  behavior 

Tectonic  strain  release  at  the  site  was  not  revealed 
by  teleseismic  data  for  LONG  SHOT  although  close-in  seismic 
measurements  indicated  some  non-circularity  in  the  energy 
radiation  at  the  source.  As  will  be  shown  in  a  later  section 
of  this  report,  teleseismic  data  from  MILROW  does  indicate 
strain  release  or  some  other  shear-generating  mechanism 
at  the  source.  Close-in  instruments  recorded  significant 
horizontally-polarized  shear-wave  motion  (Orphal  et  al.,  1970). 
McKeown  et  al.  (1970)  report  that  some  surface  fractures 
appeared  after  MILROW  and  that  vertical  displacements  up  to 
tens  of  centimeters  were  caused  by  MILROW  out  to  distances 
of  at  least  3  km.  With  the  large  number  of  near-normal 
faults  on  the  island  (Figure  1),  permanent  vertical  displace¬ 
ments  of  this  magnitude  are  not  unexpected.  McKeown  et  al., 
however,  do  not  report  horizontal  displacements,  which  are 
excellent  indicators  of  tectonic  strain  release  (Bucknam,  1969) 


LOCATION  USING  TELESEISMIC  P  ARRIVALS 


Arrival  times  for  18  Long  Range  Seismic  Measurements 
(LRSM)  stations;  for  the  center  seismometers  at  the  VELA 
arrays  LASA,  UBO  and  TFO;  and  for  many  additional  stations 
globally  distributed  were  available  to  locate  MILROW.  The 
arrival  times  of  P  waves  recorded  at  the  LRSM  stations  and 
VELA  arrays  were  all  measured  at  the  SDL.  All  of  the  P 
arrival  times  which  had  less  than  a  seven-second  residual 
( Jeff reys-Bul len  table)  were  extracted  as  listed  from 
Earthquake  Data  Report  No.  649-69  (U.S.  Department  of 
Commerce,  ESSA,  22  October  1969).  Some  additional  arrival 
times  were  determined  at  the  SDL  after  more  records  became 
available  through  the  NOS.  In  the  following  discussion, 
all  locations  were  computed  using  the  1968  P  travel-time 
tables  (H68)  as  presented  in  the  Bulletin  of  the  Seismo- 
logical  Society  of  America,  Volume  58,  No.  4.  Only  stations 
between  20°  and  100°  epicentral  distance  were  used. 

MILROW  was  first  located  without  restraining  the  depth 
using  the  158  arrival  times  listed  in  Table  II,  The  result 
was  a  23,2  km  hypocenter  shift  at  an  azimuth  of  356°  from 
the  true  location  and  a  computed  depth  of  66  km.  This  compares 
with  the  LONG  SHOT  shift  of  29.6  km  at  an  azimuth  of  6°  with 
a  78  km  depth  using  329  stations  (also  between  20°  and  100° 
distance)  with  the  H68  travel-time  table.  Azimuthal  coverage 
for  both  events  was  excellent,  except  for  the  southeast 
quadrant. 

The  MILROW  event  was  relocated  again  using  the  158  arrival 
times  but  with  depth  restrained  to  its  actual  value.  The 
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resultant  shift  from  the  actual  location  was  22.2  km  on  an 
azimuth  of  344°  as  depicted  in  Figure  2.  This  shift  is  nearly 
coincident  with  the  depth-restrained  shift  of  20.1  km  on  an 
azimuth  of  346°  for  LONG  SHOT  using  its  329  stations,  also 
depicted  in  Figure  2.  Nearly  the  identical  location  bias  has 
appeared  for  MILROW  as  for  LONG  SHOT.  A  discussion  of  this 
bias  and  analysis  of  the  LONG  SHOT  travel-time  data  can  be 
found  in  Lambert  et  al.  (1969).  Herrin  and  Sorrells  (1969), 
Davies  and  McKenzie  (1969),  and  Jacob  (1970)  have  proposed 
structural  models  below  Amchitka  Island  to  account  for  this 
bias. 

Next,  a  network  of  73  stations  which  recorded  both  LONG 
SHOT  and  MILROW  and  which  supplied  good  azimuthal  coverage  was 
used  to  locate  both  events  with  depths  restrained  to  actual 
valo'^i.  Travel-time  residuals  calculated  for  MILROW  at  these 
73  stations  are  presented  in  Figure  3  vs  epicenter-station 
azimuth;  this  shows  the  trend  of  the  residuals  with  azimuth 
which  produces  the  bias  of  the  LONG  SHOT  and  MILROW  locations. 
With  this  common  network  LONG  SHOT  shifted  22.7  km  on  an 
azimuth  of  335°  and  MILROW  shifted  21.3  km  on  an  azimuth  of 
333°  as  shown  in  Figure  2.  The  near  identity  in  shift  is  not 
surprising  since  Figure  4  shows  anomalies  (travel-time  residuals 
referenced  to  the  UBO  residual  of  -.02  seconds)  to  be  nearly 
identical  for  the  two  events;  within  reading  error,  the 
common  network  solutions  shown  in  Figure  2  exhibit  the  same 
separation  as  the  actual  locations. 

Finally,  MILROW  was  located  applying  the  LONG  SHOT  anomalies 
shown  in  Figure  4  to  the  MILROW  arrival  times,  and  vice-versa, 
in  the  manner  described  by  Chiburis  (1960)  using  his  program 
SHIFT,  again  with  depths  restrained  to  the  known  values.  MILROW 


shifted  1.2  km  on  an  azimuth  of  21°  and  LONG  SHOT  shifted 
0.9  km  on  an  azimuth  of  190°  as  depicted  in  Figure  2.  The 
two  vectors  are  very  short  and  nearly  mirror  images  of  one 
anotlier.  Thus,  we  conclude  that  the  KILROW  and  LONG  SHOT 
sites  can  be  adequately  handled  by  a  single  set  of  travel¬ 
time  anomalies  and  that  location  bias  does  not  change'  from 
one  site  to  the  other  for  the  network  used. 
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RECORDED  PHASES 


Ampl i tudes 

AH  the  amplitude  data  that  was  available  to  the  Seismic 
Data  Laboratory  is  listed  in  Tables  III  and  IV.  Corresponding 
travel  times  are  figured  for  the  first  visible  emergence  of  , 

the  phases  recorded  on  short-period  instruments,  or  long-period 
instruments  if  the  phase  is  unrecorded  on  short-period  ones. 

Table  III  contains  data  from  the  LRSM  stations  and  VELA  arrays, 
and  this  data  was  all  visually  measured  at  the  Seismic  Data 
Laboratory.  In  Table  III,  "SPR"  and  "LPR"  refer  to  true  radial 
components  for  MILROW;  likewise,  "SPT"  and  "LPT"  refer  to 
true  transverse  components.  The  horizontal  recordings  at  most, 

MILROW  stations  were  rotated  on  the  digital  computer  to 
provide  these  oriental  ons.  Table  IV  contains  data  from  sta¬ 
tions  reporting  to  the  NOS.  Some  of  this  data  relating  to 
body-wave  and  surface-wave  magnitude  was  reviewed  at  the 
Seismic  Data  Laboratory,  but  most  of  the  information  was  * 

communicated  through  the  NOS  and  has  not  been  reviewed  by 

them  nor  verified  by  us.  All  the  amplitude  data  is  plotted  in  > 

the  amplitude-distance  graphs  that  follow  and  is  employed  in 
average  magnitude  determination. 

Several  phases  were  recorded  from  MILROW  which  were  not 
seen  from  LONG  SHOT.  These  comprised  direct  shear  waves,  long- 
period  compressional  waves.  Love  waves,  and  some  late-arriving 
core  phases.  These  were  all  of  such  amplitude  that,  with  the 
yield  difference  between  MILROW  and  LONG  SHOT,  we  presume 
that  similar  phases  for  LONG  SHOT  were  not  visible  on  tele- 
se.’smic  recordings  simply  because  they  were  below  the  noisj 
level  in  amplitude.  Figures  5  through  8  are  plots  of  recorded 
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amplitudes  vs  distance  for  the  most  widely-recorded  phases 
from  MILROW;  P,  PcP,  LQ  (Love),  and  LR  (Rayleigh), 
respectively.  "Amplitude"  here  means  one-half  the  maximum 
peak-to-peak  amplitude,  corrected  for  system  response  and 
divided  by  the  period  of  the  maximum  motion.  For  body-wave 
phases,  maximum  motion  is  measured  no  farther  than  three  or 
four  cycles  after  first  motion;  and  for  surface-wave  phases, 
maximum  motion  is  measured  on  group  arrivals  having  a  oeriod 
of  about  twenty  seconds. 

No  detailed  discussion  of  compressi onal -wave  and  Rayleigh- 
wave  amplitudes  is  desirable  since  it  would  be  repetitious  with 
the  thorough  report  on  LONG  SHOT  data.  We  have  examined  the 
MILROW  data  in  regard  to  amplitude  variations  and  have  detected 
no  change  in  the  overall  pattern  of  MILROW  P,  PcP,  and  LR 
amplitudes  compared  to  the  same  LONG  SHOT  amplitudes.  For 
body  waves,  the  data  base  is  much  less  than  for  LONG  SHOT;  and 
for  Rayleigh  waves,  even  though  several  new  stations  recorded 
these  for  MILROW,  the  new  data  is  so  sparsely  distributed  or 
else  so  near  LONG  SHOT  stations  that  it  does  not  add  any 
signif'-ant  regional  information. 

Shear  energy  and  possible  mechanisms 

The  recording  of  horizontally  polarized  shear-wave  motion 
(SH)  cannot  be  accounted  for  by  modeling  the  explosion  as  a 
spherically  symmetric  compressive  source.  However,  the  presence 
of  prominent  LQ  waves  from  NTS  explosions  have  led  investi¬ 
gators  to  propose  a  mechanism  of  tectonic  strain  release  (Brune 
and  Pom'^roy,  1963;  Toksoz  et  al.,  1965).  Kim  and  Kisslinger 
(1967)  have  shown  with  model  experiments  that  SH  motion  is 
generated  in  prestressed  media  by  explosions  and  that  only  P 
motion  is  generated  when  the  media  is  not  prestressed.  Also, 
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Kisslinger  et  al,  (1961)  have  shown  that  horizontally 
polarized  shear  motion  can  occur  due  to  crack  formation 
about  an  explosive  source.  We  cannot  ascertain  the  exact 
generating  mechanism  with  the  teleseismic  data  available, 
and  our  intent  in  this  section  is  merely  to  call  attention 
to  the  presence  of  horizontally  polarized  shear-wave  energy 
on  MILROW  recordings.  Several  stations  in  Tables  III  and  IV 
recorded  LQ,  and  group  travel  times  were  such  that  these 
waves  must  have  originated  at,  or  near,  the  MILROW  site. 
Direct  shear  waves  were  visible  on  short-period  instruments 
at  KN-UT,  SHL,  COL,  TUC,  and  6UA;  and,  although  separation 
into  radial  and  transverse  components  was  impractical  for 
the  WWNSS  stations,  rotation  at  KN-UT  positively  revealed 
a  horizontally  polarized  shear  arrival.  Long-period,  direct 
shear  waves  were  found  on  FB-AK,  KN-UT,  LC-NM,  RK-ON,  and 
NP-NT  recordings;  horizontally-polarized  shear  waves  were 
positively  identified  at  three  of  these  (FB-AK,  KN-UT,  and 
LC-NM)  after  horizontal  components  were  rotated  into  radial 
and  transverse  components.  The  horizontal  shear  amplitudes 
will  be  discussed  later  in  relation  to  their  utility  for 
discrimination.  Vertically  polarized  shear  energy  was 
recorded  at  some  of  the  above  stations,  but  this  can  be 
explained  by  incidence  of  compress ional  waves  on  any 
elastic  boundary  near  the  source  rather  than  by  some  shear¬ 
generating  mechanism  at  -the  source. 
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MAGNITUDE 


Teleseismic  magnitude  estimates  of  MILROW 

Body-wave  and  surface-wave  magnitude  estimates  of  MILROW 
were  made  according  to  routine  SDL  procedure.  Body-wave 
magnitude  for  stations  from  16°  to  100°  distance  was  computed 
by  the  formula 

nib  "  ^°3l0  ^ 

where 

A  =  zero-to-peak  ground  motion  (millimicrons), 

T  =  period  (seconds), 

B  =  distance  correction  factor  of  Gutenberg  and  Richter 
(1956). 

To  obtain  A,  the  maximum  peak-to-peak  amplitude  of  the  P  phase 
within  the  first  three  or  four  cycles  on  the  record  is  halved 
and  corrected  for  system  magnification  at  the  period  T  of 
the  maximum  amplitude. 

Figure  9  shows  m^^  vs  distance  for  the  51  stations  in 
Tables  III  and  IV  which  have  P  amplitudes.  The  average  iDj^, 
obtained  by  averagi ng the  51  individual  magnitudes  without 
weights,  is  6.42;  and  the  standard  deviation  about  this 
average  is  0.40. 

Rayleigh-wave  magnitudes  for  stations  out  to  100° 
distance  were  computed  by  the  formula: 

Ms  =  log^Q  (A/T)  +  1.66  log^^A  -  0.18 
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where 


A  =  peak-to-peak  ground  motion  (millimicrons), 

T  =  period  (seconds), 

A  =  epicentral  distance  (degrees). 

The  value  A  is  obtained  by  measuring  the  maximum  peak-to-peak 
amplitude  on  the  record  at  a  period  T  of  about  twenty  seconds 
and  correcting  for  system  magnification  at  T  seconds.  Figure 
10  shows  Mg  vs  distance  for  the  47  stations  having  LR  ampli¬ 
tudes  in  Tables  III  and  IV.  The  average  Mg,  obtained  as  for 
m^j,  is  4.84;  and  the  standard  deviation  is  0.28. 

The  fact  that  the  M„  estimate  has  less  variance  than  the 

s 

m^j  estimate  was  also  true  for  LONG  SHOT.  As  for  amplitudes,  it 
would  be  repetitious  with  the  analysis  of  LONG  SHOT  magnitudes 
to  discuss  MILROW  magnitudes  vs  distance  or  region  since  we 
examined  the  data  and  saw  that  it  essentially  duplicated 
the  LONG  SHOT  data. 

Relative  amplitude  of  MILROW  and  LONG  SHOT 

For  purposes  of  determining  a  yield  difference  and  of 

perhaps  detecting  a  change  in  the  shape  of  the  source  spectra, 

it  is  important  to  determine  accurately  the  magnitude  difference, 

both  m.  and  M^ ,  between  MILROW  and  LONG  SHOT.  To  do  this  for 
b  s 

mjj,  we  measured  the  amplitude  of  the  first  compressional  move¬ 
ment  from  zero-to-peak  for  the  two  explosions  at  common  stations 
because  this  first  quarter-cycle  of  motion  should  be  uncontami- 
nateci  with  other  arrivals  (the  surface  reflection  at  the  source 
being  the  primary  problem)  and  therefore  clearly  indicative 
of  relative  source  magnitude.  Use  of  common  stations  eliminates 
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propagation  effects.  Only  film  was  employed  in  this  analysis; 
and  for  various  reasons,  mostly  excessive  amplitudes  for 
MILROW  or  low  film  amplitudes  for  LONG  SHOT,  only  nine  stations 
provided  reliable  amplitudes  for  both  events.  The  nine  stations 
and  the  corresponding  amplitude  ratios  are  listed  in  Table  V. 
The  average  of  these  ratios  is  0.154,  which  represents  an  m^^ 
difference  of  0.81.  These  nine  amplitude  ratios  form  a  sample 
large  enough  to  determine  the  mj^  difference,  especially  in 
consideration  of  the  small  standard  deviation  (0.019)  of  the 
ratios.  We  have  neglected  period  in  these  measurements 
because  there  is  no  strong  evidence  for  claiming  that  the 
MILROW  periods  for  the  first  quarter-cycle  are  different  than 
those  of  LONG  SHOT  at  these  nine  common  stations,  and  the 
inclusion  of  a  period  measurement  to  produce  A/T  would  only 
introduce  error  into  an  already  satisfactory  measure  of 
relative  amplitudes.  Also,  on  the  assumption  of  equal  periods, 
relative  ground  motion  amplitudes  were  determined  by  dividing 
film  amplitudes  by  the  magnification  at  1.0  cps  always, 
regardless  of  actual  observed  period. 


We  can  compare  this  precisely  determined  mj^  difference 

of  0.81  between  MILROW  and  LONG  SHOT  with  the  unadjusted 

mjj  difference  of  0.57  using  all  available  data  which  resulted 

in  an  m^^  of  5.85  for  LONG  SHOT  (274  stations)  and  of  6.42 

for  MILROW  as  stated  in  the  previous  section  (51  stations). 

Thus,  routine  m^^  estimation  from  short-period  recordings 

has  produced  a  0.24  error  in  the  relative  m^^  difference 

between  MILROW  and  LONG  SHOT  if  we  accept  0.81  as  the 

correct  m.  difference.  Part  of  this  error  is  due  to  the 
b 

fact  that  different  sets  of  stations  were  used  to  estimate 
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the  magnitudes  of  the  two  explosions;  and  part  is  due  to 
surface,  and  other,  reflections  which  enhance  or  diminish 
the  amplitudes  of  peaks  after  the  first  and  distort  the 
observed  periods  of  maximum  motion. 

To  determine  relative  M  for  MILROW  and  LONG  SHOT,  we 

s 

employed  the  match-filter  concept  at  all  the  common  LRSM 
stations  (except  LC-NM  which  had  a  timing  rate  error  on 
tape  for  the  LONG  SHOT  recording  period).  LONG  SHOT  locations 
WH-YK,  CR-NB,  and  PG-BC  were  changed  to  WH2YK,  CR2NB,  and 
PG2BC  for  MILROW;  but  the  distance  of  the  move  was  so  small 
in  comparison  with  Rayleigh  wave  lengths  that  we  can  consider 
these  three  pairs  as  common  stations  for  the  two  events. 
Alexander  and  Rabenstine  (1967)  have  shown  how  the  size  of 
one  event  can  be  found  relative  to  that  of  a  reference  event 
used  to  match  filter  it  by  comparing  the  amplitude  of  the 
crosscorrelation  peak  for  the  two  events  with  the  amplitude 
of  the  autocorrelation  peak  for  the  reference  event.  We  assume 
that  the  spectra  are  identical  in  shape  so  that  the  relative 
size  of  the  two  events  remains  constant  with  frequency  for  the 
band-limited  signal  on  the  long-period  LRSM  recordings. 
Signal-to-noise  ratios  for  MILROW  are  about  10  while  for  LONG 
SHOT  they  are  typically  1  to  2.  Table  VI  lists  the  stations 
used  in  the  match  filtering  and  the  amplitude  of  the  LONG  SHOT 
LR  signals  relative  to  those  of  MILROW  obtained  by  match 
filtering.  The  .average  of  the  nine  LR  ratios  in  Table  VI  is 
0.091,  which  implies  an  M^  difference  of  1.04  between  MILROW 
and  LONG  SHOT.  For  comparison,  we  have  a  M^  difference  of 
0.78  if  we  accept  the  routine  magnitude  estimate  of  4.84  for 
MILROW  as  stated  in  the  previous  section  (47  stations)  and 
of  4.06  for  LONG  SHOT  (55  stations'.  An  examination  of  the 
visually-measured  LR  amplitudes  as  given  in  Table  III  of  this 


report  for  MILROW  end  in  Table  II  of  the  LONG  SHOT  report  reveals 
that  the  ratio  of  these  amplitudes  is  about  twice  as  high, 
on  the  average,  as  the  ratios  from  the  match-filter  output 
given  in  Table  VI.  It  is  most  likely  that  noise  on  these 
nine  LONG  SHOT  recordings  accounts  for  this  doubling  of  the 
true  LR  amplitude;  and,  if  we  can  extend  this  assumption 
to  all  other  recordings  of  LONG  SHOT,  we  can  appreciate  how 
the  real  magnitude  difference  of  1.04  is  reduced  to  0.78 
using  visual  analysis  of  LONG  SHOT  since  the  logarithm  of  2 
is  0.30.  We  conclude  that  LONG  SHOT  should  really  be  about 

=  3.80  because  we  accept  the  visually-determined  estimate 
of  MILROW  to  be  4.84  as  obtained  from  all  available  data  and 
wfi  accept  the  difference  of  1.04  as  obtained  by  match 
fi 1 teri ng. 

We  have  fixed  the  m^^  difference  between  MILROW  and  LONG 
SHOT  at  0.81  and  the  difference  at  1,04;  thus  there  is  a 
twofold  increase  in  LR  amplitudes  for  MILROW  over  what  would 
be  expected  on  the  basis  of  the  mj^  difference,  and  the  sou'^ce 
spectrum  as  a  whole  has  changed  shape  with  a  slight  boost  in 
long-period  energy  for  MILROW  compared  to  LONG  SHOT,  Source 
spectral  shapes  vs  yield  shown  by  von  Seggern  and  Lambert 
(1969)  predict  such  a  result  in  the  yieius  under  consideration 
(100  to  1000  kt);  we  reproduce  their  granite  curves  in  Figure  11. 
These  curves  indicate  that  M^  (at  Ts:20  secs)  will  always  be 
directly  proportional  to  the  logarithm  of  the  yield  while  m^^ 

(at  Txl  sec)  will  level  off  at  higher  yields  after  increasing 
proportional  to  the  logarithm  of  the  yield  up  to  about  100  kt. 
This  Mg-vs-yield  behavior  is  true  for  all  geologic  medium,  but 
the  exact  relation  of  m^  to  yield  will  vary  with  medium.  If 
we  assume  that  the  predicted  yields  of  LONG  SHOT  (80  kt)  and 
MILROW  (1000  kt)  ("re  good  estimates  of  the  actual  yields,  our 
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M  differenca  (1.04)  agrees  very'  well  with  the  difference 
s  \  • 

(1.10)  of  the  logarithms  bf  the  yields.  Also,  we  see  that 
'ith8  0.81  mjj  difference  between  LONG  SHOT  and  MILROW  agrees^ 
very  well  wi^th  the  difference  of  the  logarithms  of  the 
spectral  amplitude  at  a  period  of  one  second  for  1000  and 
80  kt  shots  in  grani te,  approx'ifmately  O.SO^in  Figure  11.  Thus  ^ 
magnitude-yield  scaling  at  Amchitka  on  the  basis  of  these 
'\two  shots  follows  closely  the  predictions  for  both  short 
Ind  long  periods  given  by  von  Sejggern  and  Lambert  who  used 
the  sourcif-spectrum  scaling  theory  of  Haskell  (1967).  Itj 
should  be  mentioned  that^ Haskell  used  empirical  measurements 
from  very  small  (<5  kt)  explosions  to  formulate  theoretical 
predictions. 

I 

Haskell  neglected  depth'.of  detonation  in  his  thepry,  but 
Mueller  and  Murphy  (1971)  have  derived  the  following  long-period 
approximation  which  includes  depth: 


where  A  is  the  amplitude,  d  is  the  depth,  and  V  is  the  yield 
of  an  underground  explosion.  This  applies  only  to  events^ 
detonated  in  the  same  medium,  and  this  condition  is  not 
strictly  satisfied  by  LONG  SHOT  and  MILROW.  Using  the  given 
values  of  depth  and  yield,  the  amplitude  ratio  A^/Ag  for  LONG 
SHOT  tip  MILROW  is  .134  which  implies  a  M^  difference  of  0.87. 
This  is  somewhat  smaller  than  the  1.04  difference  determined 
by  match  filtering  above.  Thus  Amchitka  M_-vs-yield  data 

I  5 


I 


\ 

I 
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agree  better  with  Haskell's  predictions;  however,  Isince  i 
Mueller  and  Murphy  convincingly  'establish  depth  dependency 

I 

in  their  formulajtion,  the  closer' agreement  of  Haskell's 
theory  with  the  ^ata  may  be  due  to  a  fortuitous  set  o^ 
parameters  for  the  two  shot  sites  on  Amchitka.  We  a-i^e  not 
fully  satisfied  that  a  theoretical  M^-vs-yield  rel'ation 
has  been  establisihed. 

I  / 

For  the  case'  of  mj^,  Mueller  an^4lurphy  provide  a  depth- 
'  dependent  scaling  theory  which,  aljtnough  it  is  impossi^ble  to 
assign  exact  values  to  all  the  prarameters,  indicates  an  mj^ 
difference  between  LONG  SHOT  and  MILROW  which  is  less  than 
the  0.80  predicted  by  Haskell's  simple  theory.  Again  the 
parameters  at  the  two  sites  may  be  fortuitously  valued, 
such  that  Haskell's  predictions  agree  better  with  observa¬ 
tions;  and  we  do  not  accept  as  established  aCy  m^^-vs-yield 
theoretical  scaling  relation,  \  ^ 
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APPLICATION  OF  DIAGNOSTIC  CRITERIA  TO  MILROW  AND  ITS 
COLLAPSE  WITH  REFERENCE  TO  LONG  SHOT  AND  EARTHQUAKES 


The  purpose  of  this  section  is  to  compare  the  MILROW 
explosion  and  its  collapse  to  the  LONG  SHOT  explosion,  and 
to  earthquakes,  using  diagnostic  criteria  developed  primarily 
from  NTS  explosions  and  from  earthquakes  in  the  Western  United 
States.  These  diagnostics  placed  LONG  SHOT  in  the  explosion 
population  although  severe  travel-time  bias  caused  such  great 
errors  in  epicenter  position  and  depth  that  location  by  itself 
would  not  have  identified  it  as  an  explosion.  MILROW  suffered 
the  same  location  bias  (Location  Using  Teleseismic  P  Arrivals 
Section),  However,  the  surface-reflected  pP  phase  for  both 
MILROW  and  LONG  StIOT  indicate  shallow  source  depths  in  close 
agreement  with  the  known  depths  of  detonation.  These  depths 
suggest  explosion  sources.  All  the  diagnostics  used  for  LONG 
SHOT  place  MILROW  in  the  explosion  population.  The  same 
diagnostics  tend  to  place  the  MILROW  collapse  in  the  earthquake 
population.  If  one  knew  that  LONG  SHOT  was  an  explosion,  then 
MILROW  is  positively  identified  as  such  because  of  its  epi¬ 
center  and  depth  position  determined  using  LONG  SHOT  anomalies 
and  because  of  its  close  similarity  to  LONG  SHOT  in  its  entire 
seismic  signature.  We  will  now  discuss  in  detail  the  applica¬ 
tion  of  the  various  diagnostic  criteria  to  MILROW  and  its 
collapse.  Basic  seismic  data  on  the  collapse  are  given  in 
Table  VII. 

Location  and  depth  of  focus 

Since  location  and  depth  of  focus,  determined  independently 
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of  the  seismic  signature,  are  such  unique,  and  powerful  dis¬ 
criminants,  we  consider  these  first. 

Using  raw  arrival  times,  MILROW,  just  as  LONG  SHOT,  was 
located  at  sea  and  well  into  the  upper  mantle  (Locations  Using 
Teleseismic  P  Arrivals  Section);  this  result  alone  would  cause 
one  to  suspect  it  was  an  earthquake.  If  travel-time  anomalies 
determined  from  LONG  SHOT  are  employed,  the  MILROW  location  is 
only  1.2  km  from  ground  zero.  However,  even  with  the  LONG  SHOT 
anomalies,  a  depth-free  run  of  SHIFT  with  the  73  stations  used 
in  the  Locations  Using  Teleseismic  P  Arrivals  Section  resulted 
in  a  -10.6  km  depth  (minus  meaning  above  the  surface)  compared 
to  a  66  km  depth  using  raw  arrival  times.  Thus,  with  precisely- 
known  time  anomalies,  the  depth  of  MILROW  is  not  estimated  as 
accurately  as  the  areal  location,  although  the  -10.6  km  depth 
certainly  implies  an  event  close  to  the  surface.  Large  depth 
and  origin  time  instability  in  location  is  acknowledged  for 
raw  arrival  times;  for  the  case  of  MILROW  some  improvement  in 
this  error  has  resulted  from  the  application  of  LONG  SHOT 
travel-time  anomalies.  Chiburis  and  Ahner  (1970)  have  documented 
the  depth  estimate  improvement  found  at  NTS  using  relative 
travel-time  anomalies  on  many  events,  and  this  improvement  at 
Amchitka  is  about  the  same. 

For  the  collapse,  only  six  P-wave  ar'^ival  times,  these 
being  dubious,  were  available  for  location.  Even  using  LONG 
SHOT  anomalies,  the  collapse  located  77  km  r.riay  from  ground 
zero,  at  a  depth  of  86  km. 

Cohen  (1970)  has  shown  how  detonation  depth  of  a  nuclear 
explosion  can  be  determined  from  the  spectra  of  the  P  arrival 
and  P  coda.  This  method  yielded  a  depth  estimate  of  LONG  SHOT 
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only  30  per  cent  greater  than  the  actual  depth  of  702  meters. 

The  method  involves  averaging  P-wave  spectra  from  as  many 
stations  as  possible  to  suppress  noise,  and  then  using  the 
average  spectrum  to  obtain  the  delay  time  associated  with  the 
surface-reflected  pP  by  measuring  the  periodic  undulation  in 
the  spectrum  due  to  pP-P  interference.  The  output  presented 
here  consists  of  the  average  power  spectrum  from  many  stations 
(individual  spectra  are  normalized  to  1.0  before  summing)  and 
the  dot-product  of  the  pseudo-autocorrelation  and  cepstrum 
functions.  Results  for  MILROW,  MILROW  collapse,  and  LONG  SHOT 
are  shown  in  Figure  12.  The  final  spectra  were  formed  by 
averaging  17  LRSM  and  VELA  station  spectra  for  MILROW,  six 
spectra  for  MILROW  collapse,  and  27  spectra  for  LONG  SHOT. 
Further,  the  average  spectra  from  17  sub-array  center  instru¬ 
ments  at  LASA  for  MILROW,  MILROW  collapsCv  and  LONG  SHOT  are 
shown  in  Figure  13,  The  most  prominent  negative  peak  that 
appears  in  the  dot-product  trace  at  time  x  provides  an 
objective  measure  of  the  periodic  spectral  minima  due  to 
the  interference  of  pP  with  P.  This  delay  time  x  is  related 
to  depth  by 

t7 

d  =  — - 

2 

where  V"  is  the  average  compressional-wave  velocity  of  the  shot 
overburden.  For  the  LONG  SHOT  site,  ?  was  determined  to  be 
3.48  km/sec  from  an  uphole  time  of  0.202  sec  reported  by  Day 
and  Murrell  (1967).  For  the  MILROW  site,  an  uphole  time  of 
0.320  sec  (W.R.  Perret,  Sandia  Corp.,  personal  communication) 

^or  the  explosion  gives  an  average  overburden  velocity  of 

3.73  km/sec,  which  agrees  well  with  the  bore-hole  velocity  survey 
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(Snyder,  1969),  From  Figure  12,  t  =  0.65,  0,65,  and  0.50  for 
MILROW,  MILROW  collapse,  and  LONG  SHOT  respectively.  Thus  the 
depth  for  MILROW  and  MILROW  collapse  can  be  estimated  as 
1.210  km  which  is  nearly  equal  to  the  actual  depth  of  1.216  km. 
For  LONG  SHOT  the  depth  estimate  with  x  =  0.50  is  0.870  km,  24% 
greater  than  the  actual  depth  of  0.702  km.  This  latter  result 
disagrees  slightly  with  the  result  given  in  the  LONG  SHOT  report 
because  of  certain  refinements  in  the  cepstral  analysis  and 
the  fact  that  previously  x  =  0.55  secs  was  used  for  the  LONG  SHOT 
delay  time.  The  delay  time  x  can  be  estimated  to  no  greater 
than  0.05  sec  since  the  digitizing  rate  for  the  time  series 
is  20  samples  per  second,  and  in  fact  the  choice  between  0.50 
and  0.55  is  somewhat  arbitrary.  In  Figure  13,  which  shows  the 
analysis  for  LASA  sub-array  centers,  x's  of  0.60,  0.60,  and 
0.50  second  a.-e  indicated  for  MILROW,  MILROW  collapse,  LONG 
SHOT  respectively.  The  values  of  x  using  LASA  only  for  the 
three  events  as  seen  in  this  figure  are  identical  or  nearly 
identical  to  those  using  all  LRSM''and  VELA  station  (including 
LAO)  spectra  averaged  together  as  shown  in  Figure  12. 

Since  the  cepstral  analysis  leading  to  the  dot-product 
trace  is  merely  a  means  of  detecting  the  undulation  in  the 
spectra  caused  by  pP  interference,  it  is  not  necessary  to 
utilize  it  if  the  troughs,  especially  the  first  one,  are 
readily  apparent  in  the  spectra.  Cohen  emphasizes  that  where 
X  is  of  the  order  of  0,5  or  less,  the  first  spectral  trough 
will  be  at  2.0  cps  or  greater  with  further  periodic  troughs 
hidden  by  noise;  and  the  dot  product  estimate  of  x  may  be 
unreliable.  In  the  case  of  MILROW  and  LONG  SHOT  the  first 
spectral  troughs  are  readily  apparent  in  Figure  12  at  about 
1.2  cps  and  1.8  cps  respectively.  Then,  since  x  =  1/f,  delay 
times  of  .83  and  ,56  seconds  result  for  MILROW  and  LONG  SHOT 
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respectively.  These  give  depths  of  1.550  km  for  MILROW  and 
0.970  km  for  LONG  SHOT,  which  are  actually  poorer  estimates 
than  obtained  from  the  dot-product  trace.  We  have  no  good 
explanation  of  why  the  pP  analysis  has  given  depths  which 
exceed  the  true  depths  here.  For  the  collapse,  the  summed 
spectrum  in  Figure  12  shows  no  apparent  trough  and  no  energy 
beyond  1.4  cps;  this  means  the  depth  estimate  from  the 
cepstral  analysis  is  suspect. 

One  further  aspect  of  the  depth  discriminant  is  the  lack 
on  MILROW  recordings  of  observed  depth  phases  distinctly 
arriving  several  seconds  after  P  motion.  Identification  of 
such  a  phase  would  indicate  a  crustal  or  upper  mantle  earth¬ 
quake.  However,  many  shallow  earthquakes  do  not  produce  observ¬ 
able  pP  phases,  and  absence  of  pP  alone  cannot  be  used  to 
classify  an  event  as  an  explosion. 

All  location  and  depth  information  taken  together  is 

i ncofl^-liis i ve  i n . idont-i-fy-i-n^-MILROW  as  an  explosion.  Without 

the  benefit  of  LONG  SHOT  calibration,  the  MILROW  depth  estimate 
of  66  km  implies  an  earthquake  even  though  spectral  analysis 
supports  about  a  1.5  km  source  depth.  Spectral  analysis  in 
itself  cannot  be  decisive  since  occasional  earthquakes  at  all 
depths  may  have  mechanisms  that  would  produce  phases  with 
small  time  delays  relative  to  P  right  at  the  source.  But  the 
-10.6  km  depth  from  relative  travel-time  location  and  the 
success  of  spectral  analysis  in  finding  pP  corresponding  to 
a  near-surface  event  together  lend  credence  to  the  idea  that 
this  new  event  on  Amchitka  Island  was  an  explosion  if  one  had 
no  independent  information  on  MILROW. 
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Using  the  MILROW  magnitude  estimates  made  with  all  the 
available  stations  (Magnitude  Section)  of  M^  =  4.84  and 
mjj  =  6.42  and  the  MILROW  collapse  average  magnitudes  from 
Table  VII  of  M^  =  4.06  and  m^^  =  4.20,  these  M^  vs  mj^  points 
are  compared  in  Figure  14  with  those  from  LONG  SHOT  and 
shallow  Aleutian  Island  earthquakes.  Also  included  in 
Figure  14  is  a  M^  vs  m^^  relation  (Lambert,  1971)  for  39  NTS 
events  which  illustrates  the  absolute  necessity  of  regionalizing 
the  Mg  vs  mjj  discriminant.  The  important  fact  is  that,  regard¬ 
less  of  their  relation  to  NTS  shots,  Amchitka  shots  separate 
clearly  from  shallow  earthquakes  in  the  surrounding  region. 

This  result  would  not  be  changed  if  the  mj^  of  MILROW  were 
changed  to  agree  with  the  relative  mj^  differences  from  LONG 
SHOT  as  obtained  in  the  Magnitude  Section;  in  fact,  such 
correction  would  enhance  MILROW's  separation  from  the  earth¬ 
quakes  in  Figure  14.  As  stated  in  the  Magnitude  section  we 
regard  M^  of  MILROW  to  be  wel 1 -determi ned  because  of  the 
high  signal -to-noise  ratios.  In  the  same  section,  it  was 
shown  that  the  M^  of  LONG  SHOT  should  be  revised  downward; 
this  would  further  separate  LONG  SHOT  from  earthquakes. 

The  MILROW  collapse  with  an  Mg  equal  to  that  of  LONG  SHOT 
and  an  mj^  over  an  order  of  magnitude  lower  is  in  the  far 
extreme  of  the  earthquake  population  of  Figure  14.  An  expla¬ 
nation  of  the  Mg-mjj  character  of  collapses  is  the  prolonged 
mechanism  of  these  events  which  is  explained  in  detail  by 
Houser  (1969).  He  suggests  a  time  duration  of  "several  seconds 
to  several  tens  of  seconds",  depending  on  the  medium,  for  the 
process  of  collapse.  Smith  (1963)  and  Toksoz  et  al.  (1964) 
both  conclude  from  far-field  seismic  evidence  that  the  rise¬ 
time  and  duration  of  the  collapse  displacement-time  function 
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is  greater  than  that  of  explosions.  This  mechanism  would 
tend  to  produce  small  body  waves  since  short-period  energy 
is  generated  incolierently  over  several  seconds,  but  long- 
period  energy  is  essentially  in  phase  and  additive  even 
if  the  duration  is  several  seconds. 

Shear  waves 

Short-period  SH/P  ratios  (horizontally-polarized  shear 
A/T  to  compressional  A/T)  determined  for  MILROW  from  four 
stations  give  an  average  ratio  of  0.16  (Table  VIII).  Shillong, 
India,  recorded  an  SHj  but  the  P  amplitude  was  too  large  to 
be  read.  For  LONG  SHOT,  the  one  observation  of  a  short-period 
SH  phase  was  at  Shillong;  and  this  gave  a  SH/P  ratio  of  .076. 
Thus  no  direct  comparison  between  MILROW  and  LONG  SHOT  is 
possible  for  short-period  SH  excitation. 

Long-period  SH  and  SV  were  observed  at  five  stations  for 
MILROW,  but  only  three  of  these  stations  recorded  long-period 
P  motion.  For  these  three,  long-period  SV/P  ratios  range  in 
value  from  1.30  to  5.17;  and  the  single  available  long-period 
SH/P  ratio,  at  FB-AK,  is  2.79  (Table  VIII),  No  long-period 
shear  phases  were  observed  on  LONG  SHOT  recordings,  and  so 
no  comparison  can  be  made  for  long-period  SH  excitation. 

For  MILROW,  twelve  stations  from  Table  III  had  observable 
LQ  waves,  and  the  LQ/LR  average  ratio  from  these  (Table  IX) 
is  0.56.  (Note  that  no  LQ/LR  ratio  could  be  formed  at  FB-AK 
because  LR  was  clipped.)  For  LONG  SHOT  and  the  MILROW  collapse, 
Love  waves  were  detected  and  measured  by  the  match-filter 
technique.  This  will  be  discussed  more  fully  in  the  Tectonic 
Strain  Release  Section  in  conjunction  with  the  question  of 
tectonic  strain  release.  Basically,  the  results  gave  LQ/LR 


-24- 


ratios  for  LONG  SHOT  and  the  collapse  which  were  about  equal 
to  and  about  one-third  of,  respectively,  the  ratio  for 
MILROW. 

For  comparison  with  MILROW  long-period  and  short-period 
SH/P  and  LQ/LR  ratios.  Tables  VIII  and  IX  give  similar  ratios 
taken  from  stations  recording  BENHAM  and  BOXCAR,  two  large- 
yield  NTS  underground  explosions;  these  values  show  no 
significant  difference  between  the  Amchitka  and  NTS  areas 
in  regard  to  the  amount  of  horizontally-polarized  shear 
energy  generated  relative  to  compressional  energy.  Movement 
along  known  faults  has  been  triggered  by  both  these  NTS 
explosions  (Mckeown  and  Dickey,  1969)  implying  strain  release 
as  the  mechanism  for  shear  excitation.  McKeown  et  al.  (1970) 
report  smaller  observed  movements  along  surface  faults  and 
fractures  for  MILROW  than  for  NTS  events,  but  large  (tens  of 
centimeters)  elevation  changes  indicative  of  normal  faulting 
at  depth.  So  tectonic  strain  release  may  have  played  a  role 
in  shear-wave  generation  for  MILROW  also. 

It  remains  to  compare  MILROW  SH/P  and  LQ/LR  ratios  with 
earthquakes.  We  have  for  MILROW  a  short-period  SH/P  ratio  of 
0.16,  an  average  from  four  stations;  a  long-period  SH/P  ratio 
of  2.79,  from  a  single  station;  and  a  LQ/LR  ratio  of  0.56, 
from  eleven  stations.  Unpublished  data  at  the  SDL  taken  from 
LRSM  van  and  VELA  observatory  earthquake  bulletins,  com¬ 
prising  data  on  thousands  of  earthquakes,  shows  that  recordings 
at  vhese  stations  have  on  the  average  a  short-period  SH/P  of 
2.5,  a  long-period  SH/P  of  2.0,  and  a  LQ/LR  of  0.7.  Thus 
MILROW's  short-period  SH/P  ratio  is  over  an  order  of  magnitude 
smaller  than  typical  earthquakes  while  the  long-period  SH/P 
and  LQ/LR  ratios  are  not  significantly  different  from 
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corresponding  earthquake  average  ratios.  Still,  the  long- 
period  P  and  S  phases  from  MILROW  are  about  an  order  of 
magnitude  smaller  than  similar  phases  for  earthquakes  o: 
equivalent  mj^,  as  are  the  LR  waves. 

To  regionalize  the  shear  discriminant,  we  measured 
shear  and  compressional  phases  on  recordings  of  the  VELA 
observatory  UBO  from  over  a  hundred  earthquakes  having 
epicenters  in  the  Aleutian  Islands.  We  found  that  the 
average  short-period  SH/P  ratio  was  1.21  for  earthquakes, 
about  one-half  the  world-wide  average  given  above.  Compare 
this  again  with  the  average  short-period  SH/P  from  four 
stations  of  0.16  for  MILROW.  We  mention  that  no  single 
earthquake  of  those  analyzed  had  a  short-period  SH/P 
ratio  <  0.16.  Thus  short-period  SH/P  offers  a  good  diagnostic 
for  identifying  explosions  in  the  Aleutian  region.  With  the 
same  data  base  of  Aleutian  Island  earthquakes,  an  average 
long-period  5H/P  ratio  of  1.86  emerged,  about  the  same  as 
the  world-wide  average  given  above;  the  MILROW  long-period 
SH/P  ratio  of  2.79  (one  station)  is  larger  and  does  not 
provide  any  diagnostic  aid. 

Comp! exi ty 

The  ratio  of  energy  in  the  P-wave  coda  to  that  in  the 
first  few  seconds  after  the  intitial  P  impulse  was  suggested 
by  Carpenter  (1965)  as  an  aid  to  discrimination.  The  SOL 
has  routinely  calculated  this  ratio,  termed  "complexity", 
for  many  events  over  the  past  several  years;  and  a  detailed 
explanation  of  this  calculation  is  found  in  the  LONG  SHOT 
report,  MILROW  complexities  (F^)  for  seventeen  stations  are 
plotted  as  a  function  of  distance  in  Figure  15.  Complexities 
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for  LONG  SHOT  and  the  MILROW  collapse  at  stations  common  with 
any  of  these  seventeen  are  also  plotted  in  Figure  15.  All 
these  complexity  values  are  listed  in  Table  X. 

Comparison  of  ten  common  station  complexities  for 
MILROW  and  LONG  SHOT  shows  similar  complexity  values.  The 
ten-station  average  complexity  for  MILROW  is  3.34  and  for 
LONG  SHOT  is  3.16.  These  averages  include  complexity  values 
from  stations  at  slightly  different  locations:  WH2YK,  PG2BC  , 
and  CR2NB  for  MILROW  vs  WH-YK,  PG-BC, and  CR-NB  for  LONG  SHOT. 
We  cannot  account  for  the  large  difference  in  complexities 
at  BE-FL  between  MILROW  and  LONG  SHOT. 

Comparison  of  six  common-station  complexities  between 
MILROW  and  MILROW  collapse  shows  that  the  collapse  is  much 
more  complex  than  the  explosion.  The  common-station  average 
complexity  for  MILROW  is  2.55  and  for  the  collapse  is  6.25. 

Thus,  as  expected,  the  complexities  for  MILROW  are 
approximately  equivalent  to  those  of  LONG  SHOT.  Further, 
the  collapse  is  about  2.5  times  more  complex  than  either 
MILROW  and  LONG  SHOT;  this  is  explained  by  our  previously 
stated  (in  the  Application  of  Diagnostic  Criteria  Section) 
concept  of  the  collapse  mechanism  which  must  be  of  greater 
time  duration  to  account  for  much  greater  energy  arriving 
in  the  P  coda  than  for  the  explosion.  The  LONG  SHOT  average 
complexity  (2.51)  is  about  one-half  that  of  an  Andreanoff 
Island  earthquake  (4.53)  recorded  at  common  stations  as 
determined  by  Lambert  et  al.  MILROW,  with  an  average 
complexity  nearly  equal  to  LONG  SHOT,  could  be  distinguished 
from  this  earthquake  also  on  this  basis;  but  the  MILROW 
collapse  average  of  6.25  would  place  it  above  this  earthquake 
in  complexity. 


P-wave  spectra 


Comparison  of  the  P-wave  spectra  for  MILROW,  MILROW 
collapse,  and  LONG  SHOT  at  common  stations  should  exhibit 
prima^'ily  differences  in  source  functions  and  depth  since 
these  three  events  share  the  same  or  nearly  the  same  epi¬ 
center  and  thus  hava  identical  path  and  site  effects  on  their 
signals.  One  factor  in  the  scaling  and  identification  problem 
is  the  relative  content  of  longer-period  to  shorter-period 
energy  as  a  function  of  magnitude.  However,  interference  by 
pP  causing  spectral  minima,  as  discussed  in  the  Application 
of  Diagnostic  Criteria  Section,  can  complicate  the  P-wave 
spectra. 

To  obtain  spectra  for  this  report,  Fourier  amplitude 
coefficients  are  computed  for  a  signal  and  a  noise  sample  of 
equal  length;  then  the  noise  and  signal  coefficients  are 
squared  and  the  noise  spectrum  subtracted.  The  resultant 
spectrum  is  reduced  to  an  amplitude  spectrum  corrected  for 
system  response  and  static  magnification.  The  spectrum  is 
then  plotted  for  a  frequency  band  of  0,2  to  6.0  cps  in  terms 
of  acceleration  density  (my/sec  /cps).  Signal  sample  lengths 
range  from  2,5  to  5.3  secs  (real  time)  after  the  first 
detectable  motion  of  P,  Spectra  for  MILROW  and  LONG  SHOT 
for  nine  common  stations  are  shown  in  Figure  16.  Nine 
additional  spectra  for  MILROW  are  shown  in  Appendix  I. 

Six  spe.*-ra  for  the  MILROW  collapse  are  shown  in  Figure  17. 
These  figures  also  show  the  corresponding  signal  transformed. 

Comparison  of  the  waveforms  and  spectra  among  events  at 
common  stations  show  significant  differen  ‘  even  though 
path  and  recording  site  are  nearly  identi  .  These  results 
clearly  indicate  substantial  differences  in  the  source-time 
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function*  dspth  of  detonation,  or  geological  structure  at  the 
event  sites. 

Spectral  amplitude  ratios  for  MILROW  and  MILROW  collapse 
are  determined,  as  for  LONG  SHOT,  thus: 

^2  ^4 

S,  -  I  A(f.}/^  A(f.) 

^1  ^3 

A{f^)  =  value  of  the  amplitude  spectrum  at  each  0.1  cps, 

•  f^  =  0.5  cps, 

f2  =  1.0  cps, 

f^  =1,1  cps, 

f^  =  2.0  cps. 

These  ratios  are  listed  in  Table  XI  and  show  more  energy  for 
MILROW  in  the  longer-period  portion  of  the  spectra  relative 
to  LONG  SHOT.  Further,  the  collapse  ratios  show  about  the 
same  amount  of  longer-period  energy  as  MILROW.  However,  the 
signal-to-noise  ratios  for  the  collapse  P-waves  are  low,  and 
the  spectra]  estimates  are  not  as  accurate  as  those  for 
MILROW  and  LONG  SHOT.  Comparing  ratios  indicates  that  they 
are  higher  by  a  factor  of  two  to  three  for  MILROW  except  at 
PG-BC  and  KN-UT.  Here  PG-BC  for  LONG  SHOT  is  high  due  to  an 
anomalously  low  signal-to-noise  ratio.  We  cannot  explain  the 
relationship  of  ratios  at  KN-UT  other  than  to  suggest  that 
the  P  amplitude  for  MILROW  may  have  exceeded  the  linear 
portion  of  the  dynamic  range  of  the  instrument  since  the 
P-wave  magnitude  is  nearly  an  order  of  magnitude  greater 


-29- 


than  the  mean  for  MILROW  (Table  III).  It  is  clear  that  MILROW 
does  show  more  longer-period  energy  in  the  P-waves  relative  i 

to  LONG  SHOT  on  the  basis  of  the  S,  ratio  alone. 

'  \  , 

Figure  18  compares  the  average  spectra  at  18'  stations  for 
MILROW  and  27  stations  for  LONd  SHOT.  To  obtain  the  average, 
the  basic  amplitude  spectra  were  squared  and  normalized  so 
that  the  shape  of  the  spectrum  at  each  station  would  be 
given  equal  weight.  These  normalized  spectra  were  summed, 
multiplied  by  frequency  squared,  and  divided  by  the  number  of 

I 

stations  .to  obtain  thelaverage  energy  spectrum  for  each  event. 
In  Figurel  18  we  point  out  the  ^presumed  periodic  spectral 
troughs  associated  with  oP  interference;  since  data  above 
4  cps  may  merely  be  earth  noise  or  spectral  leakage,  troughs 
at  these  hi gher  \f requencies  may  be  spurious.  At  the  low 
frequency  end,  because  the  short-period  spectiral  estimates 
are  determined  from  signal  samples  of  about  3.0  to  5.0  secs, 
th^  spectra  are  generally  not  reliable  below  0.33  cps. 
Theoretical  source  energy  spectra  for  granite  and  tuff  at 
100  and  1000  kt  using  scaling  parameters  and  equations  from 
von  Seggern  and  Lambert  (1969)  are  shown  in  Figure  19,  We 
include  spectra  for  MILROW  and  LONG  SHOT  formed  by  connecting 
peaks  in  the  energy  spectra  in  Figure  18  with  smooth  lines; 
these  spectra  are  positioned  on  the  abscissa  of  the  graph 
relative  to  each  other  and  to  the  theoretical  spectra  only 
for  comparison  of  spectral  shapes.  The  peak  frequency  for  ^ 
MILROW  is  lower  than  Vhat  for  LONG  SHOT,  Frequency-dependeht 
attenuation  would  be  expected  to  shift  the  real  peaks  to 
lower  frequencies  than  the  theoretical  ones,  and  in  the  case 
of  LONG  SHOT  the  nulling  at  1.85; cps  due  to  the  interference 
of  pP  could  also  give  an  apparent  shift  to  a  lower  frequency. 
Below  1  cps  the  slopes  of  the  observed  spectra  roughly 
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correspond  to  the  predicted;  however,  for  the  higher  frequency 
portion  (1  to  4  cps)  they  do  not,  and  attenuation  increasing 
wi th frequency  should  account  for  this.  Beyond  3  cps  much  of 
the  energy  may  simply  be  leakage  f|^om  the  spectra  at  lower 
frequencies  because  of  the  rectangular  time  window  used,  iln 
ariy  case  the  peak  frequencies  do  show  the  predicted  shift  as 
a  function  of  magnitude  or  yield;  this  shift  has  been 
previously  documented  for  lower  yields  by  Berg  ]and  Papageorge 
(1964). ■  ] 

Re^garding  the  identification  problem,  short-period 
spectral  ampl'iitude  ratibs  for  many  Asian  earthquakes  and 
explosions  have  been  determined  at  LASA  by  Lac6ss  (1970) 
using  ten-second  samples.  These!  ratios  are  computed  as  follows 

rl.95  rO.85  1 


C  !  .  fU.  OO 

s,  =  A(f)df/  A(f)df 

Jl.45  JO.35 


llhe  MILROW  ratio  computed  on  this  basis  is  0.90  with  a 
magnitude  of  6.72  at  LAO,  and  the  MILROW  collapse  ratio  is 
0.16  with  a  LASA-beam  magnitude  of  4.22.  We  show  both  these 
ratios  plotted  with  those  of  Lacoss  in  Figure  20.  MILROW 
clearly  falls  into  the  explosion  (jopulation  as  did  LONG  SHOT, 

and  the  colMpse  fal'is  well  into  tihe  earthquake  poiiulation. 

'  \ 

Therefore,  MILROW  P-wave  spectra  behave  as  expected  for 
explosions.  The  collapse  ratio  at  LASA  sl^ows  the  collapse  to 
be  typical  of  earthquakes;  but  as  Table  XI  shows  the  six 
collapse  spectral  ratios  are  highly  variable,  and  the  low 
signal-to-noise  values  for  this  event  makes  spectral . rati  os 
and  thus  discrimination  based  on  them,  dubious. 
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Rayleigh-wave  spectra 

von  Seggern  and  Lambert  (1969)  have  indicated  that  for 
explosions  the  shape  of  the  theoretical  source  spectrum  in 
the  frequency  band  of  0.067  cps  (T  =  15  sec)  to  0.020  cps 
(T  =  50  sec)  should  not  change  with  magnitude.  Therefore 
MILROW  and  LONG  SHOT,  having  virtually  the  same  epicenter, 
should  have  identical  Rayleigh-wave  spectra  at  common  stations. 
With  regard  to  collapse  mechanisms,  Houser  (1969)  suggested 
a  time  duration  of  several  seconds  to  tens  of  seconds  which 
is  dependent  upon  the  medium.  There  is  evidence  from  the 
short-period  P-wave  spectra  at  LASA  and  LRSM  stations 
(Figures  12  and  13)  that  more  longer-period  energy  is 
present  for  the  collapse  relative  to  MILROW.  Further,  the 
Mg-vs-m^j  relationship  also  implies  that  the  source  function 
is  extended  over  a  greater  span  of  time  than  for  explosions. 
Thus,  we  expect  more  longer-period  energy  in  Rayleigh-wave 
spectra  for  the  MILROW  collapse  relative  to  MILROW  and  LONG 
SHOT. 

Rayleigh-wave  amplitude  spectra  are  estimated  in  the 
same  manner  as  described  in  the  preceding  section  on  P-wave 
spectra  for  velocity  windows  of  about  3.7  to  2.5  km/sec. 
Spectral  estimates  for  twelve  stations  common  to  MILROW  and 
LONG  SHOT  are  shown  in  Figure  21,  again  in  terms  of  acclera- 
tion  density.  Seven  additional  spectra  for  MILROW  are  shown 
in  Appendix  II.  The  shapes  of  the  spectra  are  very  similar 
for  five  of  the  twelve  common  stations  (PG-BC,  KN-UT,  RK-ON, 
TFO,  HN-ME)  for  MILROW  and  LONG  SHOT  in  the  frequency  band 
0.020  to  0.067  cps;  however,  as  stated  previously,  the 
s i gnal -to-noi se  ratios  for  LONG  SHOT  ranged  from  one  to  two 
and  thus  the  LONG  SHOT  spectra  are  greatly  affected  by  the 
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noise.  The  average  normalized  energy  spectra  for  the  twelve 
common  stations  in  Figure  22  show  that  MILROW  and  LONG  SHOT 
are  similar  within  the  frequency  band  of  about  0.02  to  0.075  cps. 
(Average  normalized  energy  spectra  are  determined  as  described 
in  the  preceding  section.)  At  higher  frequencies,  0.075  to 
0.10  cps,  the  noise  dominates  the  spectra  for  LONG  SHOT  and 
comparison  with  MILROW  is  impossible. 

Thirteen  Rayleigh-wave  spectral  es<’imates  for  the  MILROW 
collapse  are  shown  in  Figure  23.  Here  most  of  the  spectra 
correspond  to  MILROW  with  respect  to  the  positions  of  spectral 
maxima  and  minima,  but  most  also  show  much  more  longer-period 
energy.  The  average  normalized  energy  spectra  for  MILROW  and 
the  collapse.  Figure  24,  clearly  illustrates  this  fact. 

Energy  ratios  between  adjacent  bands  of  frequency  for 
Rayleigh  waves  have  been  used  by  Lambert  and  von  Segqern  (1969) 
as  an  aid  in  the  identification  problem.  For  events  recorded  at 
distances  greater  than  1000  km  and  M^  greater  than  3.0,  their 
preferred  ratio,  R^ ,  is  determined  as  follows: 


Rl  =  - 


-  I  (f  E„(f)df/f  E„{f)df) 
n  ^  'J.  m'  m'  ' 

m=l  '1  '2 


where  R^  equals  the  average  energy  ratio  for  n  stations  record¬ 
ing  one  event  and  is  the  energy  spectrum  at  station  m.  The 
limits  of  integration  are: 
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=  0.02C8  (T^  =  48  seconds), 
fg  =  0,0455  (Tg  =  22  seconds), 
fj  =  0.0667  (Tj  =  15  seconds). 


Calculated  ratios  for  MILROW,  LONG  SHOT,  and  the  MILROW 
collapse  are  plotted  as  a  function  of  magnitude  in  Figure  25. 
Explosion  ratios  are  independent  of  magnitude  only  if  the 
signal-to-noise  ratio  is  high  enough  to  allow  a  valid 
analysis.  The  Rayleigh-wave  magnitude  for  this  figure  is 
determined  from  the  spectral  estimates  by  having  the  computer 
pick  the  maximum  A/T  in  the  ground  velocity  spectrum  between 
periods  of  17  and  23  seconds  and  use  this  value  in  the 
standard  formula  given  in  the  Magnitude  Section.  There¬ 
fore  these  magnitudes  will  usually  differ  somewhat  from 
those  determined  from  film  analysis. 

In  Figure  25  R.j  determined  from  19  stations  for  MILROW 
is  clearly  in  the  explosion  population  which  forms  a  horizontal 
band  at  »  0.1.  In  addition,  the  value  is  nearly  equal  to 
that  of  LONG  SHOT.  IT-j  determined  from  13  stations  for  the  MILROW 
collapse  is  greater  than  for  MILROW  and,  in  fact,  greater 
than  for  any  Aleutian  earthquake  analyzed.  Thus  R^  shows 
that  the  observed  MILROW  Rayleigh-wave  spectrum  in  the  15  to 
50  seconds'  band  is  similar  to  other  explosions  while  the 
collapse  spectrum  is  definitely  not. 

Although  the  Rayleigh-waves  at  individual  stations  common 
to  both  MILROW  and  LONG  SHOT  do  not  show  exactly  identical 
spectra,  mostly  due  to  the  low  signal-to-noise  ratios  for  LONG 
SHOT,  the  average  spectral  ratios  do  show  that  the  energy  in 
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adjacent  bands  of  frequency  (0.028  to  0.0455  and  0.0455  to 
0.0667)  are  about  equal.  With  regards  to  the  collapse  and 
MILROW,  the  similarity  in  positions  of  maxima  and  minima  at 
individual  stations  plus  similarity  of  spectral  shapes  at 
shorter  periods  indicates  similar  radiation  patterns  between 
events.  The  presence  of  more  longer-period  energy  in  the 
collapse  signals  suggests  a  longer  time  function  or  larger 
spatial  dimensions  for  the  source.  We  can  certainly  eliminate 
the  second  possibility  from  mere  physical  intuition  of  the 
collapse  mechanism  and  from  its  smaller  magnitude  as  reflected 
by  or  m,^. 

Radiation  patterns 

Amplitude  and  phase  radiation  patterns  of  seismic  energy 
from  the  source  have  been  utilized  very  little  in  discrimina¬ 
tion.  One  major  difficulty  is  that,  in  determining  radiation 
patterns  for  short-period  or  long-period  waves,  the  influence 
of  path  and  recording  site  on  amplitudes  must  be  minimized.  We 
do  not  feel  that  we  have  sufficiently  accurate  information  on 
these  factors  to  undertake  the  proper  amplitude  equalization 
of  body  and  surface  waves  from  MILROW.  Figure  26  displays  all 
the  available  m^^  values  from  Tables  III  and  IV  vs  epicenter- 
station  azimuth  in  a  polar  plot;  Figure  27  does  likewise  for 
Mg.  The  distance-correction  factors  used  in  determining  m|^ 
and  Mg  according  to  the  standard  formulas  given  in  the  Magnitude 
Section  are  the  only  amplitude  equalization  applied.  It  is 
apparent  that,  in  addition  to  the  erratic  scatter  of  the  data, 
the  azimuthal  coverage  is  unsatisfactory  in  both  cases;  and 
no  'adiation  pattern  cai:  be  inferred  with  confidence.  The 
systematically  low  values  of  Mg  in  the  azimuth  range  of  60°-90° 
(western,  central,  and  southern  United  States  stations),  also 
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observed  for  LONG  SHOT,  are  just  as  probably  . the  result  of 
path  and  site  effects  as  of  a  radiation  pattern  at  the  source. 

In  the  LONG  SHOT  report  P  and  LR  radiation  patterns  of 
LONG  SHOT  were  fitted  to  theoretical  patterns  based  on 
double-couple  force  systems,  and  it  was  concluded  that  a 
circular  radiation  pattern  is  as  valid  a  fit  as  any  non¬ 
circular  pattern  associated  with  double-couple  mechanisms. 

Since  we  have  examined  the  MILROW  data  in  relation  to  LONG  SHOT 
with  polar  plots  of  amplitude  ratios  for  body  and  surface  waves 
at  all  stations  common  to  both  events  and  have  seen  no  good 
indication  of  a  different  energy  radiation  oattern  for  MILROW, 
we  will  extend  this  conclusion  to  MILROW  without  performing 
the  programmed  radiation-pattern  fitting.  We  show  in  Figure  28 
the  polar  plot  of  LR  and  LQ  ratios  determined  by  match  filtering 
the  LONG  SHOT  recordings  with  the  MILROW  signals  as  described  in 
the  Magnitude  Section  and  Application  of  Diagnostic  Criteria  to 
MILROW  and  its  Collapse  with  Reference  to  LONG  SHOT  and  Eartn- 
quakes  Section,  respectively,  and  P-wave  first-quarter-cycle 
ratios  determined  by  visual  measurement  as  described  in  the 
Magnitude  Section.  Since  only  one  quadrant  of  this  type  of 
data  is  available,  any  conclusions  are  entirely  speculative; 
however,  we  see  no  strong  indication  of  radiation  pattern 
differences  between  LONG  SHOT  and  MILROW  here. 

Other  investigators  (Brune  and  Pomeroy,  1963;  Toksoz 
et  al.,  1965)  have  used  explosion-collapse  ratios  at  NTS 
to  delineate  better  the  explosion  radiation  patterns  with 
the  assumption  that  the  collapse  released  no  tectonic  strain 
and  so  had  a  circular  radiation  pattern.  We  do  not  attempt  this 
for  MILROW  because  again  the  available  values  woulC  be  confined 
to  less  than  one  quadrant  and  also  because  we  have  already  found 
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by  match  filtering  LQ  waves  that  the  MILROW  collapse  did 
generate  horizontally-polarized  shear  waves  and  was  there¬ 
for  not  an  event  with  a  circular  pattern. 

One  further  aspect  of  radiation  patterns  is  first  motion, 
which  can  only  reliably  be  used  to  identify  earthquakes  and 
not  explosions.  The  absence  of  raref actional  first  motions  for 
a  suspicious  event  may  only  mean  that  the  recording  coverage 
was  incomplete  rather  than  that  it  was  an  explosive  source;  on 
the  other  hand,  the  presence  of  clear  raref acti onal  first 
motions  for  an  event  positively  identifies  it  as  an  earthquake. 
Evernden  (1969)  has  pointed  out  that  the  lack  of  rarefactional 
first  motions  outside  USSR  and  China  is  characteristic  of 
Kamchatka/Kuril  earthquakes  because  of  the  orientation  of  the 
fault  planes  and  the  slip  vectors.  No  rarefactional  first 
motions  were  observed  for  MILROW;  however,  for  the  reason 
just  given,  this  is  only  a  feeble  diagnostic  aid,  and  its 
improvement  as  a  diagnostic  aid  for  this  source  region  requires 
the  study  of  many  earthquakes  in  the  area  surrounding  Amchitka. 
P  waves  from  the  MILROW  collapse,  recorded  at  only  six  stations 
(Table  VII),  were  too  low  in  s i gnal -to-noi se  ratio  to  determine 
direction  of  first  motion. 

Radiation  pattern  analysis,  including  first  motion,  will 
be  a  marginal  aid  to  discrimination  in  the  Aleutian  source 
region  until  precise  information  on  path  effects  can  be 
assembled,  more  earthquake  mechanism  solutions  for  this  region 
become  available,  and,  perhaps  most  important,  better  azimuthal 
recording  coverage  with  LRSM-quality  stations  is  supplied. 

Synopsis  of  diagnostics 

We  have  applied  the  principal  diagnostic  criteria  to  the 


MILROW  event  and  to  its  collapse,  and  we  will  summarize  this 
investigation  in  a  tabular  format  where  each  diagnostic  pro¬ 
vides  a  decision  on  whether  the  event  was  an  explosion.  Most 
diagnostics  did  not  provide  unambiguous  answers  and  this  will 
be  indicated  by  question  marks.  Application  of  diagnostics  on 
a  regional  basis  is  more  meaningful  than  on  a  global  basis,  and 
we  have  attempted  to  regionalize  discrimination  for  M^-vs-mj^ 
and  for  shear  excitation  by  presenting  the  results  for  many 
Aleutian  earthquakes.  All  the  other  diagnostic  aids  suffer 
from  lack  of  Aleutian  earthquake  analysis  in  respect  to  these 
diagnostics;  however,  such  analysis  would  require  a  great 
expenditure  of  time.  For  LONG  SHOT,  one  shallow  earthquake 
of  equivalent  mj^  in  the  Andreanoff  Islands  was  studied  in 
detail  and  compared  with  LONG  SHOT;  we  will  weigh  the  answer 
to  some  of  the  classification  questions  by  comparison  of  MILROW 
and  its  collapse  to  LONG  SHOT  and  by  the  differences  between 
LONG  SHOT  and  this  one  earthquake. 

Table  XII,  then,  is  the  synopsis  of  the  diagnostic  criteria. 
It  shows  that  with  the  exception  of  location  and  depth,  MILROW 
is  identified  as  an  explosion  by  all  the  criteria  although  only 
two  are  considered  as  unambiguous.  With  the  assistance  of  LONG 
SHOT  travel-time  calibration,  the  location-depth  criteria  would 
also  tend  to  classify  it  as  an  explosion.  Also,  with  the 
knowlege  that  LONG  SHOT  was  an  explosion,  one  would  more 
positively  identify  MILROW  as  an  explosion  via  the  other 
diagnostics  because  of  the  similarity  of  these  two  events  over 
the  entire  seismic  signature  at  all  common  stations.  For  the 
collapse,  mostly  earthquake  characteristics  are  indicated  by 
the  diagnostics.  The  most  powerful  discriminant,  M^-vs-m^^,  very 
definitely  places  it  outside  the  explosion  population.  It  has 
a  complexity  2.5  times  that  of  MILROW  and  greater  than  the 
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selected  Andreanoff  Islands  earthquake  previously  mentioned. 
Short-period  and  long-period  spectral  ratios,  of  lesser  value 
In  discrimination,  place  It  In  the  earthquake  population. 
Radiation  pattern  Information  for  the  collapse  Is  uniformative 
or,  at  best,  vague.  In  other  words,  the  MILROW  collapse  would 
definitely  not  have  been  designated  an  explosion  but  may  have 
been  designated  as  an  earthquake  from  analysis  of  teleselsmic 
data  alone.  Only  the  inability  to  locate  it  accurately  In 
relation  to  MILROW  because  of  small  P-wave  amplitudes  at  tele- 
seismic  distances  prevents  Its  designation  as  a  collapse, 
barring  Independent  information  on  it,  since  as  a  whole  Its 
seismic  signal  characteristics  were  as  expected  from  a  cavity- 
collapse  mechafilsm. 
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TECTONIC  STRAIN  RELEASE  ON  AMCHITKA  ISLAND 


It  has  been  shown  that  seismicity  in  the  Aleutians  as 
seen  by  a  local  network  is  nearly  void  of  shallow  (<40  km) 
earthquakes  along  the  island  arc  (Murdock,  1969).  At  Amchitka 
Island  in  particular,  a  very  localized  network  of  land  and 
ocean-bottom  seismometers  revealed  only  infrequent  and  low- 
magnitude  shallow  shocks  and  only  very  proximate  and  weak 
seismicity  around  the  MILROW  site  after  the  detonation 
(Engdahl  and  Tarr,  1970;  Adams  et  al.,  1970).  The  MILROW 
aftershocks  were  all  confined  to  mj^  below  3.4  (personal 
communication  from  E.R.  Engdahl,  NOS).  Morris  (1970) 
states  that  faults  on  Amchitka  Island  have  been  inactive 
for  at  least  a  quarter  of  a  million  years.  Thus  there  appears 
CO  be  little  tectonic  strain  accumulation  and  release  in  the 
crust  around  the  MILROW  site.  In  contrast,  strain  accumulation 
in  the  crust  of  the  NTS  region  is  an  accepted  fact;  release 
■>f  this  strain  by  shallow  earthquakes  is  well  documented 
(Molnar  et  al.,  1969;  Slemmons  et  al.,  1965).  Also  attesting 
to  a  higher  strain  field  at  NTS  are  the  aftershocks  of  many 
NTS  explosions,  notably  BOXCAR,  which  had  at  least  seven  after¬ 
shocks  with  mjj  >  4.0  (Boucher  et  al.,  1969).  Therefore,  we 
believe  that  another  type  of  shear-generating  mechanism  such 
as  crack  formation  as  discussed  by  Kisslinger  et  al.  (1961) 
could  have  a  significant,  if  not  the  dominant,  role  in  pro¬ 
duction  of  horizontally-polarized  shear  waves  observed  on 
MILROW  recording:.  Again,  we  must  state  that  without  thorough 
near-source  seismic  measurements  and  geologic  observations  we 
cannot  precisely  define  the  mechanism  of  shear  generation. 

One  piece  of  evidence  which  clearly  demonstrates  that 
MILROW  Love  waves  cannot  be  due  entirely  to  mode  conversion 
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is  that  for  the  collapse  the  amplitude  ratio  of  Love  to 
Rayleigh  waves  was  smaller  than  that  of  MILROW,  The  reasoning 
is  that  if  conversion  were  to  account  for  the  Love  waves,  the 
relative  conversion  to  Love  waves  should  be  equal  for  both 
the  shot  and  the  collapse.  This  argument  has  been  used  pre¬ 
viously  in  the  cases  of  HAYMAKER  and  MISSISSIPPI  by  Toksoz 
et  al.  (1965).  The  evidence  is  as  follows.  Both  LONG  SHOT 
and  the  MILROW  collapse  have  an  M^  =  4.06  from  visual  analysis. 
Thus  the  Rayleigh-wave  magnitude  difference  of  1.04  between 
LONG  SHOT  and  MILROW  obtained  from  match  filtering  (Magnitude 
Section)  also  applies  to  MILROW  and  its  collapse.  Since  Love 
waves  were  not  visible  on  LONG  SHOT  and  MILROW  collapse 
recordings,  match  filtering  was  the  only  method  to  determine 
LQ/LR  ratios.  It  was  judged  that  only  four  MILROW  LRSM 
stations  had  an  LQ  wavetrain  of  sufficient  signal-to-noise 
ratio  over  a  sufficient  duration  to  use  as  a  match  filter: 
PG2BC,  RK-ON,  KN-UT,  and  SJ-TX.  KN-UT  was  eliminated  from 
the  analysis  because  of  a  tape  recording  problem  during  the 
collapse  surface-wave  arrivals.  Both  LONG  SHOT  and  MILROW 
collapse  recordings  at  the  remaining  three  stations  were 
band-passed  from  .025  to  .065  cps  over  the  expected  LQ 
arrival  window.  Horizontal  traces  for  all  recordings  were 
rotated  into  a  purely  transverse  component.  The  outputs  of 
the  match-filter  program  for  these  three  stations  are  shown  . 
for  LONG  SHOT  in  Figure  29  and  for  the  MILROW  collapse  in 
Figure  30.  Love-wave  amplitudes  of  the  two  events  relative 
to  those  of  MILROW  were  obtained  from  the  match-filter  outputs 
and  are  listed  in  Table  XIII,  From  the  figures  it  is  evident 
that  the  LONG  SHOT  match-filter  peaks  near  the  predicted  time 
on  the  last  traces  are  more  prominent  than  those  from  the 
collapse.  (Predicted  times  for  the  match-filter  peaks  in 


this  and  subsequent  figures  are  based  on  the  origin  times 
given  In  Table  Ij.)  Note  that  the  MILROW  collapse  peaks 
appear  to  be  negative  and  about  one-quarter  to  one-eighth 
cycle  later  than  If  the  source  functions  for  Love  waves 
were  180“  out-of-phase.  A  180“  phase  reversal  Is  well- 
documented  for  Rayleigh  waves  from  explosion-collapse 
pairs  (Brune  and  Pomeroy,  1963;  Smith,  1963;  Toksoz  et  al., 

1964);  but  no  explosion-collapse  phase  relation  has  been 
previously  reported  for  Love  waves,  and  In  fact  the  authors 
are  unaware  of  any  previously  reported  detection  of  Love 
waves  from  a  nuclear  cavity  collapse.  To  affirm  this  phase 
reversal  for  long  periods  between  MILROW  and  Its  collapse, 
the  collapse  Rayleigh  waves  at  the  same  three  stations  were 
match  filtered  with  the  MILROW  Rayleigh-wave  recordings.  As 
shown  In  Figure  31,  the  peaks  are  definitely  negative  and 
again,  as  for  the  Love  waves,  approximately  one-quarter  cycle 
later  than  expected.  Similar  match  filtering  of  LONG  SHOT 
Rayleigh  waves  In  Figure  32  shows  the  expected  positive 
crosscorrelatlon  peaks  at  approximatley  the  expected  time. 

From  Table  XIII  the  average  Love-wave  amplitude  ratio  between 
LONG  SHOT  and  MILROW  Is  .109  Implying  a  Love-wave  magnitude 
difference  of  0.96.  This  Is  nearly  equal  to  the  Rayleigh-wave 
magnitude  difference  of  1,04  between  LONG  SHOT  and  MILROW  and 
is  an  interesting  and  important  result.  Whatever  the  shear¬ 
generating  mechanism  of  MILROW  was,  LONG  SHOT  had  a  proportional 
amount,  almost  exactly  scaled  to  the  yield  difference  between 
the  two.  Again  from  Table. XIII ,  the  average  Love-wave  amplitude 
ratio  between  the  MILROW  collapse  and  MILROW  Itself  is  ,035 
implying  a  Love-wave  magnitude  difference  of  1.46  which  is  .42 
greater  than  the  Rayleigh-wave  magnitude  difference  between  the 
two.  And  we  regard  the  value  ,035  tc  be  In  error  on  the  high 
side  because  a  significant  portion  of  the  amplitude  of  the  match- 
filter  peaks  in  the  correlation  traces  for  the  collapse  must 
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be  noise  contributions  (more  than  for  the  LONG  SHOT  peaks 
at  least).  Thus  the  MILROW  collapse  had  an  LQ/LR  ratio  of 
at  least  three  times  lower  than  LONG  SHOT  and  MILROW  itself, 
and  we  cannot  accept  mode  conversion  to  explain  completely 
Love  waves  from  Amchitka  events  for  this  reason.  The  answer 
to  why  the  MILROW  collapse  should  have  any  Love  waves  at  all 
may  be;  1)  some  degree  of  mode  conversion  along  the  path; 

2)  non-circularity  in  the  source  space-time  function;  3)  inho¬ 
mogeneity  in  the  immediate  source  area;  4)  peculiar  geometry 
of  bedding  planes,  joints,  and  faults  around  the  MILROW  shot 
point;  5)  asymmetric  release  of  MILROW-induced  strain  in  the 
surrounding  medium.  The  last  reason  is  suggested  because 
Engdahl  and  Tarr  (1970)  report  that  the  MILROW  collapse 
brought  an  abrupt  cessation  of  MILROW  aftershock  activity 
as  though  the  main  collapse  was  the  final  event  in  the  reduc¬ 
tion  of  MILROW-induced  strain  in  the  surrouding  medium  to  a 
value  sufficiently  small  that  no  further  sudden  strain 
adjustments  would  occur.  If  the  residual  strain  field  from 
MILROW  were  asymmetric  at  the  time  of  collapse,  Love  waves 
could  be  generated.  However,  we  do  not  have  sufficient 
evidence  to  prefer  this  explanation  for  collapse  Love  waves 
over  any  of  the  others  listed  above.  Likewise,  for  MILROW 
itself  (and  LONG  SHOT  too)  we  cannot  make  definite  conclusions 
about  the  shear-generating  mechanisms.  For  the  shots  we  can 
only  state  that  mode  conversion  cannot  account  for  all  the 
Love  wave  amplitude  and  that  tectonic  strain  release  may  not 
have  been  operative  at  all  since  the  lack  of  shallow  earth¬ 
quakes  under  Amchitka  Island, of  recent  natural  fault  move¬ 
ment,  and  of  significant  post-shot  seismic  activity  indicates 
there  may  be  no  appreciable  ambient  tectonic  stresses  in  the 
upper  crust  there. 
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RELAlivE  DELAY  TIMES  FOR  MI.LROW  COLLAPSE 
AND  LONG  SHOT  SURFACE  ''WAVES 


The  differences  between  actual  and  expecte^  times  for 
the  match-filter  output  peaks  in  Figures  29  to  32  for  LONG 
SHOT  and  the  MILROW  collapse  deserve  further^  attention. 

In  Figure  30  and  31  the  negative  collapse  peaks  are 
conspicuously  later  than  expected  if  we  assume  that  the 
source  mechanism  causes  a  polarity  reversal  relative  to  j 

MILROW  but  that  its  source  time  function  is  'nearly  equivalent. 

On  the  other  hand,  the  polarity  reversal  and  the  large  time 
delay  are  both  explained  if  we  assume  that  the  Rayleigh  waves 
are  generated  by  the  impact  of  the  collapsed  cavity  material 
on  the  floor  of  .the  cavity  because  a  downward  point  force 
generates  Rayleigh  waves  180°  out  of'  phase  with  explosion¬ 
generated  waves  (Harkrider,  1964)  and  because  the  generation 
0^  high-energy  Rayleigh  wav^es  could  be  delayed  by  the  free-fall 
time  through  the  cavi^ty  of  material  connected  with  ttie  genera¬ 
tion  of  weak  initial  P-waves  which',  had  emergent  arrivals  at 
teleseismic  distances.  The  average  delay  of  the  negative  peaks 
in  Figures  3]o  and  31  is  5.0  1,4  sec  (95%  confidence  interva'Js) 

and  the  cavity  diameter  predicted  by  the  formula  of  Closmann 
(1969)  would  be  largq  enough  to  allow  a  delay  of  this  much. 

This  same  phenomenon  of  Rayleigh-wave  delay  and  p'plari'i.y 
reversal  has  been  found  for  the  BENHAM  cavity  collapse  (von 
Seggern,  1971).  ' 

The  L0|NG  SHOT  match-filter  peaks  show  advances  from  the 
predicted  t'imes  in  Figures  29  and  32,  but  only  by  a  second  or 
two.  The  predicted  times  do  take  into  account  the  sm^ll 
differences  in  epicentral  distance  to  the  three  stations  PG-BC, 
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RK-ON,  and  SJ-TX  from  the  two  detona|tion  points.  The  average 

time  advance  for  the  LONG  SHOT  surface  waves  in  Figures  29  and 
\  32  is  1,6  +  1.0  (95%  cjonfidence  intervals).  Since  we  do  not 

doubt  the  reported  origin  times,  we  believe  some  difference 
in  the  source  functions  related  to  yield,  to  interference  of 
the  pp| phase,  or  to  contrasts  in  the  mediums  immediately  | 
surrounding  the  two  test 'si tes  lis  responsible  for  this  result. 
The  dependence  of  the  phase  spectrum  of  the  source  on  yield 
can  be  derived  from  Haskell's  (1967)  representation  of  the 
source  potential.  The  derivation  |s  straightforward,  but 
long,  and  we'  give  only  the  final  expression  for  the  phase  of 
the  Fourier  transform  of  the  displacement-time  function  at 
the  source:  '  '  i  ■ 

ia)k[(24B+l  )(5u)'’-10a)^k^+k^)-.(w^-10a)^k^+5k^)  J 
\  e(a))  =  arctan - ^ - ; - 

M^(24B+l)(u^-l,0u^k^+5k^)+k^(5uf-10r/k^+k'') 

1 

1  -1/3 

The  parameter  k  scales  as  (yield)  '  and  the  dimensionless 
'  parameter  B  i£  independent  of  yield;!  both  are  dependent  on 
medium.  We  ar^  interested  in  evaluating  the  phase  lag  at  a 
period  of  20  ^ecs  for  yields  0{f  80  and  1000  kt'in  hard  rock. 

We  used  Haskell's  values  of  k  and  B  for  granite  and  scaled  k 
from  his  values fW  5  kt.  The  phase  lag  of  LONG  SHOT  was 
calculated  to  be  .0454  radians  while  that  of  MILROW  was  .0959 
radians  at  a  period  of  20  secs;  t,his  gives  approximately  a 

!  j 

0.2  second  advanife  of  the  LONG  SHOt  waves  relative  to  those 
of  MILROW.i  This  is  then  insufficient  to  expl'iain  the  observed 
early  arrival  of  LONG  SHOT  surface  waves  relative  to  those]  of 
]  MILROW.  Thus,  Haskell's  scaling  theory  may  not  be  right  for 


A  P 


the  phase  spectrum  or  else  medium  contrasts  and  differences 
In  source  depth  may  be  causing  the  observed  phenomenon. 
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s 

«» 


-46- 


SUMMARY  AND  CONCLUSIONS 


The  proximity  of  MILROW  and. LONG  SHOT,  their  similar 
detonation  mediums,  and  their  large  yield  difference  enabled 
us  to  obtain  a  limited  magnitude-yield  relation  for  the 
Amchitka  test  area.  Surface-wave  magnitude  scaled  directly 
proportional  to  the  logarithm  of  the  announced  yields,  a 
result  predicted  by  tha  work  of  Haskell  (1967).  The  m^j 
difference,  not  a  simple  function  of  yield,  also  agreed  with 
the  predictions  of  Haskell.  Because  Haskell  did  not  include 
depth  as  a  variable  in  his  formuation,  we  regard  as  somewhat 
fortuitous  the  excellent  agreement  between  M^  and  mj^  observa¬ 
tions  and  predictions  based  on  his  theory.  The  Amchitka  scaling 
appears  to  be  parallel  to  NTS  high-yield  scaling  for  events  in 
tuff,  but  slightly  higher  mj^  and  much  lower  M^  for  a  given  yield 
on  Amchitka  are  due  to  entirely  different  travel  paths  and 
different  geological  environments  in  the  Amchitka  and  Nevada 
test  sites.  The  fact  that  dominant  periods  for  NTS  Rayleigh 
waves  in  North  America  are  from  10  to  16  sec  while  those  for 
MILROW  were  from  17  ^5  sec  jreatly  accentuates  M^  for  NTS 

events  relative  to  MILROW  as  shown  by  Evernden  and  Filson  (1971), 

Oepth-'jnrestrained  locations  of  MILROW  with  raw  arrival 
times  rev'.'aled  the  same  large  (s'ZO  km)  northward  location  bias 
and  the  same  depth  error  (s70  km)  that  attended  LONG  SHOT. 
Relative  travel-time  anomalies  were  determined  for  73  LONG  SHOT 
stations  which  also  recorded  MILROW;  application  of  these 
anomalies  reduced  the  location  error  of  MILROW  to  1.2  km  along 
the  surface  and  to  -11.8  km  in  depth  (the  location  having  been 
made  above  the  earth's  surface).  So,  although  relative  travel¬ 
time  anomalies  provide  horizontal  location  accuracy,  they  still 
fail  to  provide  a  precise  depth  estimate.  This,  however,  is  not 
surprising  because  of  the  inherent  depth  and  origin  time 
instability  in  the  location  algorithm. 
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Relative  magnitudes  for  MILROW  and  LONG  SHOT  were  found 
by  measuring  the  amplitude  of  the  first  quarter-cycle  at 
selected  common  stations  in  the  case  of  P  waves  and  by 
measuring  the  match-filter  crosscorrelation  peak  of  LONG 
SHOT  and  MILROW  and  the  autocorrelation  peak  of  MILROW  in 
the  case  of  Rayleigh  waves.  In  this  way  the  m^^  of  MILROW  was 
found  to  be  0,81  higher  than  that  of  LONG  SHOT,  a  figure 
which  is  0.24  greater  than  the  m^^  difference  found  by  averaging 
for  each  event  all  the  available  m|^  values  determined  in  the 
routine  manner.  Likewise,  the  precision  M^  difference  from 
match-filter  processing  was  1.04,  a  figure  0.31  greater  than 
the  Mg  difference  found  from  routine  visual  analysis. 

Both  short-period  and  long-period  instruments  recorded 
direct  shear  phases  from  MILROW;  horizontally-polarized  shear 
motion  was  positively  identified  in  some  cases.  Also,  MILROW 
Love  waves  were  detected  at  many  stations.  The  detection  of 
previously  unseen  Love  waves  from  LONG  SHOT  followed  when  MILROW 
recordings  were  used  as  match  filters.  It  was  surprising,  how¬ 
ever,  to  detect  Love  waves  from  the  MILROW  collapse  when  the 
same  match-filters  were  applied  because  detection  of  Love 
waves  from  collapses  has  never  been  reported.  Release  of 
MILROW-induced  strain  is  a  plausible  explanation  for  the 
collapse  Love  waves.  Also  crack  formation;  mode  conversion; 
source  asymmetries;  local  geologic  inhomogeneities;  and  the 
geometry  of  bedding  planes,  faults,  and  joints  may  account 
for  these  recorded  Love  waves  from  the  collapse.  For  MILROW 
itself,  due  to  the  absence  of  shallow  earthquakes  in  the 
Aleutian  region  and  the  limited  aftershock  activity,  we 
conclude  that  the  above  various  processes  suggested  as 
possibly  operative  for  the  collapse  caused  the  shear  motions 
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observed  and  not  tectonic  strain  release  commonly  associated 
with  NTS  shots. 

In  the  identification  problem,  MILROW  behaved  as  expected 
from  the  authors'*  previous  knowledge  of  explosion  signals  and 
in  particular  from  our  study  of  LONG  SHOT  characteristics.  We 
would  identify  MILROW  as  an  explosion,  especially  if  com¬ 
parisons  with  LONG  SHOT  are  maclo.  On  the  other  hand,  the  MILROW 
collapse  definitely  does  not  classify  as  an  explosion.  An 
adequate  summary  of  the  individual  diagnostics  has  already 
been  given  in  Table  XII. 

The  major  conclusions  of  this  analysis  are  as  follows: 

1)  The  Amchitka  source  region  has  been  adequately  cali¬ 
brated  with  respect  to  P  travel  times, 

2)  Source-spectra  vs  yield  scaling  has  followed  Haskell's 
predictions  closely.  In  determining  this  scaling,  precise 
relative  magnitude  differences  should  be  obtained  by  measuring 
initial  P-wave  excursions  at  common  stations  and  from  Rayleigh- 
wave  match-filter  outputs. 

3)  While  tectonic  strain  release  has  been  the  preferred 
explanation  for  Love  waves  from  NTS  explosions,  it  probably  does 
not  have  a  large  role  in  Amchitka  explosions.  Thus  the  shear¬ 
generating  mechanism  depends  on  the  test  area,  and  probably  no 
general  model  is  appropriate. 

4)  Diagnostic  criteria,  especially  M^  vs  m^^,  developed 
primarily  from  NTS  shots  and  Western  United  States  earthquakes 
and  from  Asian  events  have  worked  well  in  the  Amchitka  region; 
however, if  hypothetically  we  were  required  to  identify 
further  explosions  in  the  Aleutians,  a  more  comprehensive 
study  of  earthquakes  in  this  region  would  be  helpful. 
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Figure  1.  Geology  of  southeastern  part  of  Amchitka  Island. 
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Figure  10.  Surface-wave  magnitude  (M  )  vs  distance  for  MILROW. 
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Figure  12.  Depth-phase  (pP)  analysis  using  averaged  normalized 
for  MILROW  (17  stations),  MIIROW  collapse  (6  stations),  and 
LONG  SHOT  (27  stations). 
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Figure  13.  Depth-phase  (pP)  analysis  using  LASA  subarrays' 
average  spectra  for  MILROW,  MILROW  collapse,  and  LONG  SHOT. 
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Figure  15.  Complexity  values  for  MILROW,  LONG  SHOT,  and  the 
MILROW  collapse  vs  distance. 
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Figure  16.  P~wave  spectra  for  stations  common  to  MILROW  and 
LONG  SHOT. 


mjj/SEC  X  10“  mp/SEC  «  lo 


Figure  16.  (Cont'd.) 
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Figure  17.  P-wave  spectra  for  the  MILROW  collapse 
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Figure  17,  (Cont‘d.) 

K-wavs  spectra  for  the  MILROW  collapse 
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Figure  19.  MILROW  and  LONG  SHOT  smoothed  average  P-wave  spectra 
compared  to  theoretical  source  spectra. 
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Figure  20.  LASA  spectral  ratios  of  MILROW  and  its  collapse 
compared  to  data  of  Lacoss  (1970)  on  Asian  earthquakes  and 
presumed  explosions. 
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Figure  21.  (Cont'd.) 

Rayleigh-wave  spectra  for  stations  common  to  MILROW  and 
LONG  SHOT. 
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Figure  21.  (Cont'd.) 

Rayleigh-wave  spectra  for  stations  common  to  MILROW  and 
LONG  SHOT. 
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Rayleigh-wave  spectra  for  stations  common  to  MILROW  and 
LONG  SHOT. 
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Figure  22.  The  average  Rayleigh-wave  normalized  energy  spectra 
for  MILROW  and  LONG  SHOT  using  twelve  stations  common  to  both 
events. 
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Figure  23.  (Cont'd.) 

Rayleigh-wave  spectra  for  the  MILROW  collapse. 
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Rayleigh-wave  spectra  for  the  MILROW  collapse. 
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Rayleigh-wave  spectra  for  the  MILROW  coll 
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Figure  25.  Raylei gn-wave  spectral  ratio  vs  Rayleigh-wave  magnitude 
for  many  events. 
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Figure  28.  Polar  plot  of  LONG  SHOT/MILROW  LQ  and  LR  matched- 
filter  amplitude  ratios  and  first  quarter-cycle  P-wave  amplitude 
ratios  vs  epicenter-station  azimuth. 


Figure  31.  Match  filtering  of  MILROW  collapse  LR  waves  (middle 
trace)  using  MILROW  recordings  (top  trace). 
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TABLE  V 


LONG  SHOT  -  MILROW  P  Amplitude 
Ratio  From  First  Quarter-Cycles  of  Motion 


STATION 

LONG  SHOT 
MILROW 

TUC 

.179 

ALQ 

.149 

KN-UT 

.132 

HN-ME 

.lo4 

LC-NM 

.149 

RK-ON 

.167 

TFO 

.179 

UBO 

.126 

LAO 

.145 

Average  ratio  =  .154 
Standard  deviation  =»  ,019 


’  18 


TABLE  VI 

LONG  SHOT  -  MILROW  LR  Amplitude 
Ratio  From  Match-Filter  Output 


LONG  SHOT 


STATION 

MILROW 

WH-YK* 

.060 

NP-NT 

.090 

PG-BC* 

.096 

KN-UT 

.086 

RK-ON 

.118 

CR-NB* 

.116 

SO-TX 

.085 

HN-ME 

.078 

BE-FL 

.088 

Average  ratio  =  .091 
Standard  deviation  =  .014 


*Moved  to  WH2YK,  PG2BC,  and  CR2NB  for  MILROW 
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TABLE  VIII 


Summary  of  Shear-to-Compresslonal 
Ratios  for  MILROW,  BENHAM,  and  BOXCAR 


MILROW  Short-Period  Shear 

A/T 


Station 

Phase 

(mu/sec) 

S/P  Ratio 

SHL 

“''h 

53. 

P  unreadable 

COI. 

S 

223. 

.27 

QUA 

'h 

235. 

.25 

TUC 

Sh 

42. 

.06 

KN-UT 

Sh 

167. 

.05 

.16  Average 

MILROW 

Long-Period 

LC-NM 

Sh 

12.6 

No  P 

LC-NM 

Sy 

8.5 

No  P 

KN-UT 

Sh 

66.7 

No  P 

KN-UT 

Sy 

54.2 

No  P 

FB-AK 

Sh 

81.0 

2.79 

FB-AK 

Sy 

150.0 

5.17 

RK-ON 

Sy 

10.1 

1.30 

NP-NT 

Sy 

34.5 

1.40 

BENHAM 

Short-Period 

NP-NT 

Sh 

38.0 

.10 

BENHAM 

Long-Perl od 

PG2BC 

Sh 

114. 

1.48 

RK-ON 

Sh 

175. 

2.92 

FB-AK 

^h 

114. 

1.43 

BOXCAR 

Short-Period 

RK-ON 

Sh 

274. 

.12 

BOXCAR 

Long-Period 

RK-ON 

8K 

1.42 

WHO 

Sh 

125. 

1.25 

'21 


TAfiLE  n 

Summary  of  i.ove-tc-Rayleigh  Ratios 
for  MILROW,  3ENHAM,  and  BOXCAR 


MILROW 

Station 

LQ  A/T 
(mu/sec) 

Lq/LR 

WH2YK 

82.0 

.23 

NP-NT 

19.5 

.10 

PG2BC 

46.3 

.51 

LAO 

54.1 

.72 

UBO 

17.6 

.43 

KN-UT 

44.2 

.87 

RK-ON 

35.0 

.27 

LC-NM 

26.5 

.70 

BY-IO 

40.9 

.94 

WQ-IL 

25.3 

.26 

SJ-TX 

48.8 

1.14 

.56 

BENHAM 

KN-UT 

19,100 

1.30 

UBO 

4,740 

1.07 

PG2BC 

1 ,950 

.40 

BOXCAR 

AT-NV 

14,000 

.66 

EY-NV 

11,300 

.29 

BF-CL 

3,850 

.59 

WW-UT 

9,860 

.86 

KG-AZ 

5,240 

.41 

ND-CL 

25,800 

.49 

CP-CL 

13,700 

.36 

TFO 

3,520 

.20 

UBO 

3,300 

.92 

PG-BC 

2,390 

.24 

SV3QB 

713 

.46 

.59 

Average-HILROW 


Average>BENHAM  and  BOXCAR 


TABLE  X 

Complexities  for  MILROW,  LONG  SHOT 
and  the  MILROW  Collapse 


Station 

Distance 

(Deg.) 

MILROW 

LONG  SHOT 

MILROW  Collapse 

FB-AK 

21.5 

4.19 

6.78 

WH2YK 

26.7 

9.63 

9.92 

PG2BC 

34.6 

6.79 

7.93 

LAO 

47.2 

2.18 

2.14 

3.64 

UBO 

49.0 

1.68 

5.01 

KN-UT 

49.1 

2.87 

2.63 

6.00 

RK-ON 

51.5 

1.20 

1.22 

TFO 

51.8 

2.76 

2.28 

8.92 

LC-NM 

56.0 

1.61 

1.41 

7.13 

CR2NB 

56.2 

2.31 

1.66 

BY-IC 

58.6 

1.63 

WQ-IL 

61.3 

1.08 

SJ-TX 

64.4 

1.70 

1.32 

AS-PA 

66.1 

2.11 

EU2AL 

66.8 

2.66 

PJ-PA 

67.9 

2,13 

BE-FL 

73.1 

2.32 

1.10 

;  Q 


r- 


m. 


fl 


6-| 


TABLE  XI 

Short-Period  P-wave  Spectral  Ratios  for 
MILROW,  LONG  SHOT  and  the  MILROW  Collapse 


SPECTRAL  RATIO  (S^) 


Stati on 

MILROW 

LONG  SHOT 

MILROW  Co 

FB-AK 

1.39 

1.83 

WH2YK* 

2.0A 

0.84 

PG2BC* 

6.77 

11.52 

LAO 

2.16 

1.64 

2.66 

UBO 

11.37 

10.44 

KN-UT 

0.86 

2.99 

10.01 

RK-ON 

1,30 

0.52 

TFO 

2.67 

1.61 

LC-NM 

5.40 

2.67 

5.09 

CR2NB* 

1.62 

0.36 

BY-IO 

2.28 

---- 

WQ-IL 

1.19 

SJ-TX 

6.12 

1.71 

____ 

GH-MS 

5,17 

____ 

AS-PA 

1.81 

EU2AL 

2.14 

PJ-PA 

3.32 

BE-FL 

4.18 

1.84 

*Were  at 

WH-YK,  PG-BC,  . 

u  CR-NB  for  LONG  SHOT 

»  ‘  I  ^ 


Synopsis 

Diagnostic 

Location  and  Depth 
of  Focus 


Ms  "’b 


Shear  Haves 


Complexity 


Short-Period 

Spectra 


Long-Period 

Spectra 


TASLE  Xil 

of  Diagnostics  for  HILROH  and  the  HILROH  Collapse 

Does  It  Identify  Event  as  an  Explosion 

HILROH _ Collapse _  _ Comments 


No 


Yes 


Yes 


Yes? 


No  Cepstral  pP  depth  analysis  gives  good 

results  for  both.  Ultfi  LONG  SHOT 
travel-time  anomalies,  HILROH  would 
appear  as  a  near-surface  event.  Even 
with  anomalies,  collapse  locates  off 
Amchitka  Island. 

No  HILROH  Is  clearly  In  explosion  popula¬ 

tion,  and  collapse  Is  clearly  in 
earthquake  population. 

Yes?  HILROH  short-period  SH/P  ratio  Is  less 

than  any  Aleutian  earthquakes  analyzecL 
anc'  long-period  shear  waves  (direct  S 
and  LQ)  are  both  low  relative  tc  mp. 

No  visible  direct  shear  phases  from 
the  collapse  were  found. 

No?  Complexities  for  only  one  earthquake 

In  the  Aleutians  have  been  reported 
(Lambert  et  a1.,  1969). 


Yes?  No?  Only  Aslan  presumed-exploslon  and 

earthquake  ratios  were  available  as 
the  two  background  populations 
(Lacoss,  1969).  HILROH  spectra 
differed  considerably  from  those  of 
LONG  SHOT  and  Individual  spectral 
ratios  were  generally  higher  than 
for  LONG  SHOT 


Yes?  No?  HILROH  follows  trend  of  explosions 

from  three  test  areas,  but  Is  amidst 
six  Aleutian  Islands  earthquakes. 
HILROH  spectra  were  very  much  like 
LONG  SHOT.  Collapse  had  more  longer- 
period  energy. 


No  rarefactlonal  first  motion  for 
HILROH.  Analysis  of  more  Aleutian 
earthquakes  Is  necessary.  HILROH 
patterns  do  not  appear  to  differ 
from  those  of  LONG  SHOT. 


Radiation  Patterns, 
Including  First  Notion 


Yes? 


? 


1 


\ 

t 

) 


station 

PG2BC 

RK-ON 

SJ-TX 
I  AVERAGE 


\ 


TABLE  XIII 

Relative  Values  of  LONG  SHOT  and  MILROW 
Collapse  LQ  Amplitudes  to  Those  of  MIL^ROW 


LOING  SHOT  LQ 
MILROW  LQ 


.099 

.143 


.109 


MILROW  Collapse  L( 
MiLROW  LQ 


•.029 

.039 


.038 


-.035 


\ 


\ 


\ 
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MILROW 
PJ-PA 
7SS0  Km 


APPENDIX  II 

ADDITIONAL  MILROW  RAYLEIGH-WAVE  SPECTRA 
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